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INTRODUCTION 


Under contract with the Department of Energy, the National Aeronautics 
and Space Administration Is managing 48 operational solar heating and cooling 
test sites (OTS) and numerous large commercia] systems. The sites are expected 
to provide information which will help engineers and designers evaluate the 
solar system's effectiveness in (l)reducing fossil fuel consumption and 
(2)maintaining design performance goals. 

The information from these operational test sites (OTS) is provided as 
a large mass of data (insolation, flow rates, temperatures, etc.) stored on tape. 

In order to attain adequate engineering assessments of these systems, procedures 
must be developed to analyze the data on a real time basis. This report presents 
a detailed discussion of three computer programs (TAPFIL, CHPLOT, WRTCNL) developed 
by MSFC to read the OTS site data tapes and put the data in a form usable with other 
solar heating and cooling system analysis computer programs. These programs utilize 
the MSFC Systems Analysis and Integration Laboratory's PDP 11/70 computer facility. 
TAPFIL verifies the analytical methodology, CHPLOT plots the site data for examination 
of system/subsystem trends, and WRTCNL lists the data in tabular form. 

Implementation of the program code on another computer system would be subject 
to an understanding of the program execution and the software capabilities of that 
system. The computer code development procedures are presented and discussed in 
such a manner as to prepare the reader to utilize this program code on any computer 
which has adequate capabilities. 


T . DATA SOURCE 


Each OTS solar system built under MSFC management utilizes a site data 
acquisition system (SDAf?) developed by IBM. The SDAS's are designed to retrieve 
raw data from 48 site measurements every r > minutes and 20 seconds continuously. 
Several large OTS systems require more than 48 measurements of data to determine 
total system performance and require two SDAS's. This* allows a maximum of 96 
channels of site data to be retrieved. When two SDAS's are used at one site, 
two methods are used to join them together. One method is to connect the two of 
them to one clock so that the measurements may be sent to the central data 
processing facility in time sequence. The other method is to connect each SDAS 
to a separate clock, thus creating two time profiles which must be sorted and 
merged at the central data processing facility to obtain a time sequence. 

The distinction between which type of time profile is created at a site 
is made by the system subtype. There are five subtypes. The subtype is determined 
from reading the IBM 360 data tapes. These tapes contain EBCDIC characters in 
a certain order according to subtype. These characters are sequentially stored 
in 62 records (210 bytes long) per data block of 13020 bytes. Each record contains 


I 


time information, site identification information and 48 channels of data. 

The order of the data arrangement and time profile is defined below according 
to subtype. The arrangement of the data on the tape is shown In Figure 1. 


SUBTYPE 

DATA ACQUISITION 

CHARACTER REPRESENTATION 

1 

One 48 channel SDAS 

• Byte 1: No significant use 



• Bytes 2-4: 3 digit Site ID number 



• Bytes 5-16: Year, Month, Day, Hour, 
Minute, Seconds 

For example - YYMMDDHHNNSS is represented 

by Bytes 5, 6, ...15,16. If date and 

time were 9/1/79 and 10:32:58, the tape would 
yield the following Byte representation - 
790901103258 (Bytes 5-16) 



• Byte 17: System Subtype 



• Byte 18: No significant use 



• Bytes 19-210: 48 4-byte hexadecimal 

EBCDIC words representing 
channels 1-48 

(There is 

no subtype 2) 


3 

1st 48 channel 
SDAS of sorted 
and merged time 
profile 

Same as above except always used in conjunction 
with sybtype 4. Bytes 19-210 represent 
channels 1-48. 

4 

2nd 48 channel 
SDAS of sorted 
and merged time 
profile 

Same as above except always used with subtype 3. 
Bytes 19-210 represent channels 49-96. 

5 

1st 48 channel 
SDAS of connected 
time profile. 

Same as above except always used with subtype 6. 
Bytes 19-210 represent channels 1-48. 

6 

2nd 48 channel 
SDAS of connected 
time profile. 

Same as above except used with subtype 5, 
Bytes 19-210 represent channels 49-96. 
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These records are stored in a data block. Subtype 1 records, utilizing 
one SDAS to retrieve data, are stored one record at a time with each time point 
being 5 minutes and 20 seconds apart. Subtype 3 records should be read immediately 
preceedlng subtype 4 records. These two subtypes should have the exact same time 
point. The subtype 4 record should then be followed by another subtype 3 record 
with a new time point 5 minutes and 20 seconds later and continue in this subtype 3 - 
subtype 4 sequential manner until completion of 270 output records. In some Instances 
however, contiguous records of subtype 3 or subtype 4 may be encountered. In these 
Instances, the contiguous subtype records are not stored and the number of output 
records will be less than 270. Subtypes 5 and 6 are read and stored exactly the same 
as subtypes 3 and 4 except no contiguous records will be encountered. 

All the above mentioned subtypes are written on 9-track data tapes by an 
IBM 360 computer at a density of 1600 bits per inch (BPI) with one end-of-file (EOF) 
written after each month of site data. 

Each OTS has an Instrumentation Program and Components List (IPCL) which is 
used to correlate the channel number to the actual measurement at the site. Listed 
in the IPCL are the channel numbers, the measurement numbers, the operating ranges, 
and the names of each measurement. Once the IBM data tape has been successfully 
read and the data stored, the IPCL is used to determine which measurements are 
needed for calculations and plots. A sample IPCL can be found in Figure 2. The 
data read from the IBM tape is stored by channel number in a data file. This will 
allow computer programs to be developed that do certain required calculations 
dependent upon a set of input channel numbers. Since every operational test site 
has a different channel orientation, those programs can be used to analyze and/or 
evaluate OTS systems according to the input channel numbers. This type of data 
storage (direct access unformatted) was chosen because it requires very little 
disk storage space compared to other types of storage, and is easily accessible, 
in that no data formatting is required. 

We have discussed the formulation of the data source (IBM 360 OTS Tapes), 
the tape format, the storage of the source data on the PDP 11/70 system disk, 
and the use of the IPCL to determine channel information for analysis and 
plotting. With this background, the specific programs used to read, convert 
and analyze the source data will be discussed next. 



FIGURE 1. IBM 360 TAPE DATA ARRANGEMENT 
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FIGURE 2 
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XI. PROGRAM USAGE 


The following three programs are designed to be run on the MSFC PDF 11/70, 
utilizing FORTRAN IV PLUS. CHPLOT and WRTCNL can be run only after a direct access 
unformatted file has been created by TAPFIL. 

Appendix A contains a sample problem utilizing all three programs. An OTS 
site contracted by WORMSER and located in Columbia. S.C. is analyzed in this problem. 
Appendix B contains flowcharts of the three programs, utilizing only the main 
subroutines. The functions performed by these subroutines are also listed in this 
appendix. A FORTRAN progrtm listing for TAPFIL can be found in Appendix C, for 
CHPLOT in Appendix D, and for WRTCNL in Appendix E. 

A brief description concerning the operation of these programs will now be 
presented. 


A. TAPFIL OPERATIONAL PROCEDURES 


The program TAPFIL is designed to create a direct access unformatted file 
from an IBM 360 tape. This file can then be used with the channel plotting program 
CHPLOT and the tabulation program WRTCNL, using the PDP 11/70. Due to the software 
used in the program, TAPFIL must be run from a TEKTRONIX 4014 Demand Terminal. 

After program execution has begun, the first statement of Figure 3 will 
appear on the screen. The first user input will be the number of channels at the 
site. To ensure program execution, 48 is the miminum number that can be entered. 

The user should enter 48 even if the site actually has fewer than 48 channels. 

If the site has more than 48 channels, the input should be in Increments of 16, 
ranging from 48 to 96. The maximum number that can be input is 96, as this is the 
number of channels on two SDAS's. 

The second user input is the desired date of data to be retrieved. This 
should be entered as MO DA YR. The program searches through the time profile on 
the tape, finds the date Input by the user, and sets flags to begin output to the 
system disk. 

The next user input is the data file name. This name should have a maximum 
of 29 characters. The convention used at MSFC is to enter the two digit site number 
with an "s" to indicate site, then the date as MODAYR, followed by .DAT to indicate 
a data file. For example, if the user wants to use the data for June 13, 1980 from 
OTS site #12, the input will be 12S061380.DAT. 

The fourth user input will be the file number of the data date. The months 
contained on a tape are listed on the front of that tape. To obtain the file number, 
count the months from the start of the list until the desired month is reached. 

This is the file number. 


The last user Input will be a decision as to whether or not a line printer 
time profile is wanted. After all inputs are complete, an output statement will 
appear stating that the data is being retrieved. 

The next output on the screen will be a time profile, as shown in Figure 4. 

If the user has so requested, this time profile will also be output to the line 
printer. Also written will be a statement telling how many records the data is 
stored in. If an error is encountered during execution, various error messages 
will be printed to inform the user. 

s 

Once the time profile has been written, the program is in pause mode. The 
user must depreso any key to continue program execution. 

A statement will then appear requiring user input as to whether more data is 
to be retrieved. If so, the statements of Figure 3 will be repeated, with the 
exception of the headings and the question concerning the number of channels at the 
site. If no more data is to be retrieved, a statement will appear stating that the 
tape is being rewound. 

The TAPFIL program has now created direct access unformatted files. These 
files are ready for use with CHPLOT and WRTCNL to analyze the site data. 


TAPE TO FILE PROGRAM 

SYSTEMS REQUIREMENTS BRANCH 

SYSTEMS ANALYSIS AND INTEGRATION LABORATORY 
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49 

ENTER DESIRED DATE OF DATA TO IE RETRIEVED AS I NO DA VR 
M IS 79 

ENTER DATA FILE NANE 
S3SSSIS79.DAT 

ENTER THE FILE NUMBER OF THE DATA DATE . 

IS 

DO VOU WANT A LINE PRINTER TINE PROFILE 7 <1-YES,2»N0> 

DATA IS NOW BEING RETRIEVED 


FIGURE 3. TAPFIL SCREEN DISPIAY 
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B. CHPLOT OPERATIONAL PROCEDURES 


The purpose of CHPLOT is to plot the data stored in « fits created by TAPFIL. 
The program has the capability to plot as many as live plots on one graph. CHPLOT 
must be run from a TEKTRONIX 4014 Demand Terminal due to the software used In the 
program. 

When program execution has begun, the first lines of Figure 5 will appear on 
the screen. The first user input Is the file name in which the data to be plotted 
Is stored. This should be the name of a fil'- created by TAPFIL. The second user 
input is the number of plots desired on the graph. The maximum number allowable 
is five. The next two user Inputs are the captions to be written at the bottom and 
on the left margin of the plot. The first caption allows 49 characters to be output; 
the second caption allows 16 characters. 

The fifth user input is the channel numbers to be plotted, along with appro- 
priate units. The IPCL should be referenced to obtain the correct channel number of 
the measurement desired. After entering a channel number, the next input should be 
the mnemonic name and the units of that measurement. The carriage return should be 
depressed after both the name and the units. Tills is repeated until the desired 
number of channel numbers and units has been input. 

The last input shown on Fifcu.a 5 is the number of channels of data at that 
site. This is subject to the same restrictions specified in the description of 
TAPFIL. After all user inputs are completed, an output statement will appear 
stating that the data is being retrieved. 

After the data has been retrieved, a plot similar to the one in Figure 6 will 
appear on the screen. The grid is designed to fit onto a regular 8-1/2 by 11 inch 
sheet after xeroxing for use in reports to which these data are pertinent. 

When the plot is completed, the program is in pause mode. The user must 
depress any key to continue execution. The next output to the screen is shown in 
Figure 7. This display gives the user a choice of five options for plotting. 

With the exception of option 2, the scaling option, the options are self-explanatory. 
The user must make his decision, enter the number of that option, and depress the 
carriage return. 

An example of the scaling options display is shown in Figure 8. This option 
gives the user the capability of windowing in on a certain section of the plot that 
was just made. The first input to the scaling option (XMAX) is the maximum abscissa 
value to be plotted; the second is the smallest abscissa value to be plotted. The 
third and fourth Inputs are the maximum and minimum values for the ordinate to be 
plotted. After entering the last value in the scaling option, the plot in Figure 9 
will appear on the screen. The plot will have the user input maximum and minimum 
values at the appropriate locations on the grid. After this plot has been completed, 
the program is in pause mode. Depressing a key will return control to the plotting 
options in Figure 7. 
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EL 55 CHANNEL PLOTTING PROGRAM 

ENTER PLOT rilE NAME 
93S061S79.DAT 

ENTER NO. OF PLOTS ON THIS GRAPH <PAX«5) 

2 

ENTER PLOT CAPTION (49 CHARACTERS) 

UORHSER E/15/79 

ENTER CAPTION FOR LEFT SIDE OF PLOT 
TERPERATURE 

ENTER CHANNEL NO. AND UNITS 
12 T1D1 
92 T151 

ENTER NO. OF CHANNELS AT SITE 

48 

DATA is KING RETRIEVED 


FIGURE 5. CHPLOT SCREEN DTSPIAY 
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FIGURE 7. PLOTTING OPTIONS 
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FIGURE 8. SCALING OPTION DISPLAY 
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FIGURE 9. SCALED PLOT 
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C. WRTCNL OPERATIONAL PROCEDURES 


T e purpose of WRTCNL Is to list in a tabular form the data stored in a 
direct access unformatted file created by TAPFIL. The program has the capability 
to list the data for up to ten channels. Since WRTCNL does not contain any special 
routines, this program may be run from any terminal with the necessary capabilities. 

The user inputs to this program are shown in Figure 10. The first input is 
the number of channels at the site. Tills is subject to the same restrictions 
specified in the description of TAPFIL. The next input is the number of channels 
to be output, with a maximum of ten. The third input is the channel numbers to be 
output. The IPCL should be referenced to obtain the correct channel numbers of the 
measurements desired, This Is also the case for the next input, which is the measure- 
ment type and number. The last user input is the data file name of the stored data. 
This should be the name of a file created by TAPFIL. 

After all inputs are completed, a listing of the requested data will be output 
to the line printer. 


ENTER NO. OF CHANNELS AT SITE 

-48 

ENTER THE NO. OF CHANNELS TO IE OUTPUT US NAN) 

2 

ENTER THOSE CHANNELS NOS. ASl 01,02,03, 

12 92 

ENTER NEASURENENT TYPE , AND NUNDER, ASlTTlQO.EPIOI.UUIOD, 

NOTE l PROGRAN EXPECTS 2 LETTER IDENTIFICATION OF NEASURENENT TVPE- 
SXNGLE LETTER NEASURENENT TYPES SUCH AS TEMPERATURE 
SHOULD IE INPUT AS TT,UU,ETC. 

TT101,TTIS1 

ENTER THE DATA FILE NAME OF THE STORED DATA 
•3S061S79.DAT 


FIGURE 10. Wimm SCREEN DISPLAY 


III. PROGRAM SOFTWARE IMPLEMENTATION 


The previously described programs were written for use on a PDP 11/70 system 
with the 11-M operating version. Software capabilities of other computer systems 
could differ greatly; therefore, the software routines used in these programs will 
be explained In detail. 

Three types of software routines m o u .id by these programs - TEKPLOT Library 
routines, UTILITIES Library routines and IV. m Conversion routines. The first and 
largest set of routines used are the TEKPLOI routines. These routines are standard 
library software routines for the PDP 11/70. Libraries on other computer systems 
contain routines with different names that perform the same functions. Before 
these programs can be used on another computer, the TEKPLOT routines must be 
replaced by similar routines compatible to the computer used. The TEKPLOT routines 
are listed and described in Table 1. 

The second type of software routines am in the UTILITIES Library. There 
are four of these routines; three of these are tape handling routines and the fourth 
is a data conversion routine. Table 2 contains a listing and description of these 
routines . 

Bytes 19-210 of a record are converted from EBCDIC floating point words to 
PDP 11/70 words by a subroutine not included in the UTILITIES Library. This routine 
is known as CFP360 and is used as CFP360(DAT(1) ,48,0) . The first argument represents 
a real data array dimensioned large enough to contain 48 words. DAT(l) represents 
the starting position to perform the required conversion and the second argument repre- 
sents the number of words to be converted. One byte of data on the IBM 360 tape is 
equivalent to one byte on the PDP 11/70 and is represented by a L0GICAL*1 statement. 

On the PDP 11/70, one floating point word contains four bytes; therefore, 192 bytes 
of data are stored in the L0GICAL*1 data arvay from bytes 19-210. This indicates 
that 48 words are to be converted by the CFP360 routine. The last argument repre- 
sents a byte-swap operation; a byte-swap is required for conversion from EBCDIC to 
PDP words, and is accomplished by setting the fourth argument to zero. The corre- 
lation between the L0GICAL*1 data array and the real data array used in CFP360 is 
done by an EQUIVALENCE statement. These statements are used: 

DIMENSION DAT(48) 

L0GICAL*1 DATA(210) 

EQUIVALENCE (DATA( 19 ) , DAT( 1) ) 

This will equivalence the real dimensioned array DAT to the last 192 bytes of the 
L0GICAL*1 array DATA so that CFP360 can convert real data words. 

The tape handling routines are standard software routines and can easily be 
compared to routines on another computer with a tape drive, but the two data 
conversion routines, EBCASC and CFP360, were developed for the MSFC PDP 11/70. 
Therefore, the comparable conversion routine to be used on another computer would 
need to convert EBCDIC characters to compatible characters for that computer system. 
Once the use of the above routines is understood, any experienced programmer should 
be able to implement these programs on any computer using FORTRAN IV PLUS. 


14 


TABLE I 


TEKPLOT ROUTINES 

ALPHA- 

Places terminal in alphanumeric mode, cursor at current (x,y) position. 

AXIS ( XLOW , YLOW , XT JIG , YLNC , XTIC , YTTG , MARKX, MARRY )- 

Draws X and Y axes. Each line must pass through the origin. If not, 
the routine returns with no action. After the axes are drawn, the user 
is left in vector mode positioned at the origin. 

XLOW- starting position on X axis in scaled, rotated units. Must 
be the lowest point through which the X uxis will pass. 

YLOW- equivalent of XLOW pertaining to the Y axis. 

XLNG- length of X axis in scaled units 

YLNG- length of Y axis in scaled units 

XTIC- distance between tic marks on X axis. If XTIC is negative, then 
is the distance between decades for a log axis. 

YTIC- equivalent of XTIC pertaining to Y axis 

MARKX- ^O, no axis annotation for X axis 

.GT. 0 , tic mark annotation for X axis 
.LT. 0 , grid annotation for X axis 

MARRY- equivalent of MARKX pertaining to the Y axis. 

BEGIN(JBAUD)- 

Initlalizes the terminal status area and the software character generator. 
Should always be the first subroutine called in any program. 

JBAUD- the baud rate of the current installation. 


CHOUT(J)- 

Outputs a single character. 

J- ASCII value of a character to be output to the T4014 scope 

ERASE- 

Erases the screen, places the terminal in alphanumeric mode, positions 
cursor to home position. 


HOME- | 

Positions cursor to the upper left corner of the screen j 

'■i 


POINT- 


Places terminal in point plot mode and positions cursor to the current 
graphics position. 





SCALE ( XFACT , Y PACT , XORC. ,YORG) - 

Establishes X and Y scaling factors and X and Y origin 
XFACT- X scaling factor (*tekpo hits /scaled unit) 
YFACT- Y scaling factor (“tekpotnts/seaied unit) 
XORG- X origin in tckpoints 
YORG- Y origin in tckpoints 


TP AUS E- 

Causes a user to become I/O bound for an indefinite period of time. 
Tills allows the user to sit and wait while inspecting a plot without 
incurring CPU charges when working in a time sharing environment. 


TPLOT ( X , Y , I PEN , MARK ) - 

Plots an (x,y) coordinate in the current mode. 

X- X value in scaled, rotated points 

Y- Y value in scaled rotated points 

IPEN- **0, dark plot (pen tip) 

•NE, 0, bright plot (pen down) 

MARK- the number of a data mark from the following list: 

0- no data mark 

1- small x 

2- small down arrow 

3- small up arrow 

4- small square 

5- small triangle 

6- small as t crick 

11- large x 

12- large down arrow 

13- large up arrow 

14- large square 

15- large triangle 

16- large as te rick 


VECTOR- 

Draws a smooth line between the current graphics position and the (x,y) 
point specified as a parameter to any of the plotting routines. 


1 <» 






TABLE 2 


UTILITIES ROUTINES 


EBCASC ( DATA , NC) - 

Converts EBCDIC characters to ASCII characters for use on the PDP 11/70. 

DATA- name of array that contains EBCDIC characters 
NC- number of characters to be converted 

SFUN(LUN , ISTAT ,NF , IRC , IDIR) - 

Performs special functions such as rewinding, forward spacing and backspacing 
records and f i les . 

LUN- logical unit number of the tape drive 

ISTAT- status code for error interpretation; use ISTAT(l) for starting 
location 

NF- number of files on the IBM 360 tape to space 

IRC- *0, normal operation 
*1, end of file 
«2, I/O error 
*~1, end of tape 

IDIR- represents which tape function to perform 
«3, rewind 

*10, forward space files 
*13, backspace files 

TGET (LUN , ISTAT , DATA , NW) - 

Retrieves a block of data off the tape. One block contains 13020 bytes, 
or 6510 PDP 11/70 words since there are 2 bytes of L0GICAL*1 data per 
word on the PDP 11/70. 

LUN- same as above 

ISTAT- same as above 

DATA- an array data is stored in; use DATA(l) for starting location 
NW- number of words of data to be retrieved 

TWAIT ( LUN , ISTAT , NWORDS , IRC ) - 

Suspends task execution until the last I/O request for the specified LUN 
has been completed. 

LUN- same as above 

ISTAT- same as above 

NWORDS- number of words of data remaining in the block after TGET 
retrieves the required number of words 

IRC- same as above 
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Of the three programs previously described, the understanding of TAPFIL 
is the most essential. CHPLOT and WRTCNL cannot be run until TAPFIL has 
created a direct access unformatted file from the data tape. Several other 
programs which require the use of TAPFIL to retrieve the data are currently in 
use at HSFC. These programs perform various analyses pertaining to collector 
efficiency. 

CHPLOT can be used to plot dally data, thus enabling the user to find a 
day which appears to have usable data. WRTCNL can be used if a listing of this 
data is desired. After a good day has been found using CHPLOT, the data can 
then be used in the analysis of collector efficiency. 

Another program currently in the development stage at MSFC will allow 
monthly data to be plotted without the prior use of TAPFIL. This program will 
read the data straight from the IBM tape and plot it one week at a time. The 
advantage of this program will be that the user can find good days without 
using TAPFIL to create numerous individual files. After the good days have 
been located, TAPFIL can then be used to create only the necessary files. 
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APPENDIX A 
SAMPLE PROBLEM 
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TAPE TO FILE PROGRAM 

SYSTEMS REQUIREMENTS BRANCH 

SYSTEMS ANALYSIS AND INTEGRATION LABORATORY 
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TPFIL TIRE PROFILE 
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EL 55 CHANNEL PLOTTING PROGRAN 


ENTER PLOT FILE NONE 
93S9S1579.DAT 

ENTER NO. OF PLOTS ON THIS CJMPH (MAX'S) 
2 

ENTER PLOT CAPTION (49 CHARACTERS) 
WORKER 6/1S/79 

ENTER CAPTION FOR LEFT SIDE OF PLOT 
TEMPERATURE 

ENTER CHANNEL NO. AND UNITS 
12 T101 

02 T1S1 

ENTER NO. OF CHAMCLS AT SITE 

49 

RATA IS lEINO RETRIEVED 








{ 


\ 



? 


I 

f 


p 
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5 


$ 1 1 
mi 
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WRTCNL PROGRAM 


ENTER NO. OF CHAIMELS AT SITE 

41 

ENTER THE NO. OF CHANNELS TO K OUTPUT (li IWO 
2 

ENTER THOSE CHANNELS NOS. ASt 01,02.03. 

12 02 

EWER NEASURENENT TYPE . AMD NUMBER, AS:TT100,EP101,UU100. 

NOTES PR OCRAN EXPECTS 2 LETTER IDENTIFICATION OF NEASURENENT TYPE- 
SINGLE LETTER NEASURENENT TYPES SUCH AS TEMPERATURE 
SHOULD IE INPUT AS TT.UU.ETC. 

TTlOl .rriSt 

ENTER THE DATA FILE NAME OF THE STORED DATA 
03S061579.DAT 


to 

Ol 


03S061b79.DAT 




T lMt 

1 T101 

T1131 

0.0611 

99.00 

112.70 

0.1300 

99. 2o 

112.70 

0.23 o9 

93 . 8o 

112. 70 

0.3278 

93.60 

112.70 

0.9lo7 

93.66 

112.90 

0.5036 

93.9/ 

112.90 

0.5999 

93.2/ 

112.70 

kJ.68^3 

92.8/ 

112.98 

0 -/ 72 ? 

92.26 

113.30 

ID. 6611 

90.10 

113.10 

0.9500 

69.31 

113.10 

1.03d9 

09.S1 

112.90 

1.12/8 

69.31 

112.90 

1.2167 

89.31 

113.10 

1 .3036 

69.51 

112.90 

1.39*49 

69.51 

112. /B 

1.9833 

69.31 

112.30 

1.5722 

68.91 

112.50 

1.6611 

09.31 

112.30 

1 . 7500 

68.91 

112.50 

1.6389 

68.52 

112.70 

1.92/8 

68.91 

112.50 

1.9999 

08.32 

1 12.30 

2.0863 

o7.93 

112.10 

2.17/2 

67. 7^ 

112.30 

2 .26o 1 

67.73 

1 12.10 

2.3530 

67.53 

1 12.30 

2.9939 

67.33 

112.10 

2.o269 

67.33 

111.51 

2.7199 

87.73 

111.71 

2.6063 

67.33 

111.31 

2.89/2 

67.13 

111.71 

2 . 98o 1 

87.13 

111.71 

3.0730 

86.99 

111.31 

3. 1639 

87.13 

111.31 

3.2b28 

Ob . 79 

111.31 

3.6917 

8 6 . / m 

111.11 

3.9306 

66. 9*4 

110.91 

3 • b 1 99 

66.3m 

111.31 

3.o0o3 

6b . 39 

110.71 

3.69/2 

65.93 

110.91 

3.78ol 

65.53 

111.11 

3.6730 

65.73 

111.71 

3.9639 

69.70 

111.51 

9 . 0528 

62.39 

111.90 

9.1917 

bl .01 

111.71 

9.2306 

81.01 

111.31 

9.3199 

6 1 . 60 

110.91 

9.9063 

61.60 

110.91 

9.9972 

61.80 

110.91 


This is a sample of the printout produced by WRTCNL . 
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FLOWCHARTS ami REFERENCE TABLES 


WRTCNL FLOWCHART 
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CHPLOT 



CHPLOT 


WATNXT 


PSCALE 







TAPFIL SUBROUTINE LIST 


CFP360- 

CHSIZE- 

INPUT- 

RWN'D- 

TAPFWD- 

TEKPLOT- 

TIMCVT- 

UTILITIES- 


AUTOSC- 

CHSIZE- 

DATIN- 

DISPLY- 

PLOT5- 

PSCALE- 

RETREV- 

SSPLOT- 

TEKPLOT- 

INPUT- 

OUTPUT- 

TEKPLOT-- 


converts EBCDIC floating point words to PDP 11/70 words 
changes character size on a TEKTRONIX 4014 display scope 
reads user Inputs regarding site data lmformatlon 
rewinds tape 

moves tape to correct file 

routines associated with use of the TEKTRONIX 4014 scope 
converts ASCII characters to decimal equivalents 
tape bundling and data conversion routines 


CHPLOT SUBROUTINE LIST 

sets grid for the plot to the proper scale 

changes character size on TEKTRONIX 4014 scope 

retrieves floating point data from the scope 

reads user Input regarding site data information for plotting 

builds plots of data 

allows the user to plot portions of data for a particular day 
retrieves the file to be plotted from the tape 
draws axis and grid, plots and labels data 
routines associated with use of the TEKTRONIX 4014 scope 

WRTCNL SUBROUTINE LIST 

reads user input regarding site data information for printing 
prints site data information 

routines associated with use of the TEKTRONIX 4014 scope 


SO 


APPENDIX c 


TAPPTI, SOURCE LISTING 


2 THIS ROUTINE CONUERTS IBM 360 32-BIT FLOATING POINT NUNBERS 
I TO PDP-11 FORMAT. THE FORTRAN <F4P) CALLING SEQUENCE IS, 


CALL CFP360C ARRAY, N, ISUAP ) 


; WHERE, 

l ARRAY - STARTING ADDRESS OF AN ARRAY CONTAINING 

, THE 360 FORMAT DATA* IT DOESN'T HATTER IF 

, IT IS INTEGER OR REAL. IT ONLY REQUIRES THAT THE 

• 360 DATA IS IN 32-BIT PAIRS. 

I N - THE INTEGER NUMBER OF 32-BIT 360 WORDS TO CONUERT. 

I I SWAP - 0 FOR BYTE SWAP FOR DEC 

i 1 FOR NO BYTE SWAP 


.TITLE 

CFP360 


.GLOBL 

CFP360 


CFP360: 



MOW 

R0,-(SP) 


now 

R1,-(SP> 


MOW 

R2,-(SP> 


HOU 

R3,-(SP) 


HOU 

R4,-(SP) 


MOV 

RS,-(SP) 


MOO 

«4(RS),R4 

; FETCH NO. OF 32-BIT WORDS TO CONUERT 

HOU 

2(RS),R0 

; FETCH ARG. ADDRESS TO R0 

HOU 

06(R5),SWAP 

; FETCH SWAP ARGUMENT 

Sts 

CLR 

FLAG 

; CLEAR NEG. FLAG 

now 

(R0)+,R1 

{FIRST 16 BITS TO R1 

MOU 

(R0),R2 

{LAST 16 BITS TO R2 

TST 

R1 

{TEST FOR ZERO 

BNE 

?• 

; BRANCH IF NOT ZERO 


TST 

BEQ 


7f: 

TST 

BNE 

SUAB 

SUAB 

6f: 

TST 

BGE 

BIS 

If: 

nou 

CLRB 

ASL 

BIC 

BIT 

BEQ 

BIS 

2t: 

BIC 

3f s TSTB 
BLT 
ASL 
ASL 
ADC 
SUB 
BR 

4f: 

BIC 

ADD 

BIS 

8f: 

nou 

nou 

ADD 

DEC 

BGT 


R2 jTEST FOR ZERO 

8f ; BRANCH IF ZERO 


SUAP 

6 * 

R1 

R2 


; TEST FOR SUAP 
; BRANCH IF NO BYTE SUAP 
{ BYTE SUAP FOR DEC 
;SAHE 


R1 

If 

♦ 100000 , flag 


jTEST FOR NEGATIUE 
;NO 

{YES, SET FLAG 


R1,R3 

R3 

R3 

#140030, R3 
#40033,31 

2f 

♦40090, R3 


; CLEAR MANTISSA PORTION 
; SHIFT LEFT 1 
; CLEAR BIAS & SIGN 
;TEST BIAS 3IT 
jOFF 

{ON , SET IT 


#177400, R1 
R1 

4f 

R1 

R2 

R1 

#200, R3 

3f 


; CLEAR EXPONENT PORTION 
{SEE IF BINARY NORMALIZED 
; YES 

•NO, SHIFT LEFT 
;R2, SHIFT, R2 
;ADD CARRY TO R1 
; ADD BIT TO EXPONENT 


#200, R1 ;HIDDEN BIT 

R3,R1 ; COMBINE EXP. AND MANTISSA 

FLAG,R1 ;SET NEGATIUE BIT 


R2, (R0) {STORE BACK IN ARGUMENT 

R1,-(R0) 

#4,R0 ; POSITION PONTER TO NEXT 32 BIT UORD 

*4 {DECREMENT, ZERO FLAG SET IF THAT IS ALL 

SB {NO, THAT IS NOT ALL 
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cttt 

cttt 

cttt 

cttt 

cttt 


SUBROUTINE CHSIZE(N) 

THIS ROUTINE CHANGES CHARACTER SIZE ON A TEKTRONIX 4014-1 
SCOPE ACCORDING TO N AS FOLLOWS: N»1 - LARGEST 

N-2 - NEXT SNALLER 
N*3 - NEXT SNALLER 
N-4 - SMALLEST 


NP-N+55 

CALL CHOUT ( 27 ) 
CALL CHOUT CNP) 
RETURN 
END 


DISPLAY 


PIP> 


SUBROUTINE INPUT (NCH,NFILE, IIP, I REP ) 

LOG I C ALT 1 IFILE(30),IDATE(6) 

COHHON/FILE/IFILE 
CONHON/D ATE/ I DATE 
CONNON/UHCRE/XS , VS 
YS-730. 

XS-SO. 

CALL ALPHA 
LUN-7 

IF ( IREP.EQ. 1 ) GO TO 22 
CALL CHSIZE(2) 

CALL TPLOTCXS,YS,0,0) 

WRITE C I, S) 

5 FORNATdH*, 'TAPE TO FILE PROGRAN' ) 

YS-YS-25. 

CALL TPLOTCXS,YS,0,0) 

WRITE(1,6) 

6 FORNAT ( 1H* , ' SYSTENS REGUIRERENTS BRANCH') 

YS-YS-25. 

CALL TPLOT(XS,YS,0,0) 

WRITE (1,8) 

8 FORNAT (1H+, 'SYSTENS ANALYSIS AND INTEGRATION LABORATORY') 
YS-YS-30. 

CALL CHSIZEC3) 

CALL TPLOT(XS,YS,0,0) 

WRITE(1,1S) 

15 FORNATdH*, 'ENTER THE NUNBER OF CHANNELS AT SITE') 
YS-YS-15. 

CALL TPLOT(XS, VS,0,0) 

READ(1,20) NCH 
20 F0RHATCI4 ) 

VS-YS-15. 

22 CALL ?PLOT(XS,YS,t,0) 

WRITE (1,25) 

25 FORNATClN-f, 'ENTER DESIRED DATE OF DATA TO BE RETREIUED 
t AS* NO DA VR ' ) 

VS-YS-15. 


CALI TPLOT (XS,YS, 0,0) 

READ(1,30) (IDATEC I ), 1-3,6), (IDATEtl ), 1-1,2) 

30 F0RHAT(3(2A1, IX ) ) 

YS-YS-15. 

CALL TPLOT(XS,YS,0,0> 

WRITE (1# 35) 

35 FORHATdH*-, 'ENTER DATA FILE HARE " ) 

YS-YS-15. 

CALL TPL0T(XS, YS,0,0) 

READ(1,40) (IFILEd )»I>1,29) 

40 F0RHATC29A1 ) 

YS-YS-15. 

CALL TPLOT(XS,YS,0,0) 

WRITE (1,45) 

45 FORHATdH*, "ENTER THE FILE MURDER OF THE DATA DATE") 
YS-YS-15. 

CALL TPLOT(XS,YS,3,0) 

READ (1,20) NFILE 
YS-YS— 15 • 

CALL TPLOT ( XS, YS, 0, 0 ) 

WRITE (1,50) 

50 FORHATdH*, "DO YOU WANT A LINE PRINTER TIRE PROFILE 
* ? (1-YES.2-N0)" ) 

YS-YS-15. 

CALL TPLOT(XS,YS,0,0) 

RE AD (1,20) ILP 
YS-YS-15. 

CALL TPLOT (XS,YS,0,0) 

URITEd»S5) 

55 F0RHAT(1H+, "DATA IS NOU BEING RETREIUED" ) 

RETURN 

END 


PIP> 


rj </> 
or,p- 
oc ui ■ 
MCZ 

r>*-o 

mou 


jjha 

<c<rujz 

OOKlii 


38 


TAPFIL: 


C tttttttttt SECTION A **************************** 

C ********** SET UP STORAGE AND INPUTS 

DINErtSION DATH48), DAT2( 48 ) „ ISTAT ( 4 ) 

LOGICAL* 1 DATA ( 13020 ) , DATAR i (210), DATAR2 (210), DATAR ( 210 ) 
LOGICAL* 1 IFILE ( 30 ) , IDATE ( 6 ) , ICHK ( 2 ) 

COHNON/TIN/DATAR 
COHHON/FILE/IFILE 
COHNON/DATE/ I DATE 
CONNON/UHERE/XS , YS 
EQUIVALENCE ( DAT 1 ( 1 ) , DATAR 1(19)) 

EQUIVALENCE ( DAT2 ( 1 ),DATAR2( 19 ) ) 

DATA IFILE (30 )/0/ 

CALL BEGIN (9600) 

LUN-7 

CALL ASSIGN(LUN,'HT0:') 

NFILE-1 
IREP-2 
2 NRDAT*1 
- JRK-0 
YK-40. 

XH»50. 

NORECS-O 
CALL ALPHA 
CALL ERASE 

CALL SCALE(1. 0,1. 0,0. 0,0.0) 

CALL INPUT(NCH,NFILE,IL?,IREP) 

C********** ROLL TAPE FORWARD TO APPROPRIATE FILE 
CALL TAPFUD(NFILE,WFILE,ISTAT(1 ),IREP) 

NFILE-NFILE 
CALL CHSIZE(2) 

CALL ERASE 
VS-730. 

CALL TPLOT ( 305 • # YS# 0, 0 ) 

URITE(1,S) 

5 FORHAT(lH+ # 'TPFIL TINE PROFILE' ) 

CALL CHSIZE(3) 


YS-YS-1S. 

YSAV-YS 

NT-NCH+1 

cimmni open print file and storage file 

C tttttttttt ALSO PRINT LINE PRINTER HEADINGS 

IFCILP. EQ.l)OPEN(UNIT-3,NAf!E« 'PRINT. LST', TYPE- 'NEU',ERR-608) 
OPEN ( UNIT-2, NAHE- IFILE, TYPE- ' NEW ' , ACCESS- 'DIRECT ' , 
fRECORDSIZE-NT, ASSOCIATEUARIABLE-NRDAT,HAXREC-270,ERR-608 ) 
IFCILP. EG. 1) URITEC3,?) 

7 F0RHAT(1H1,T2, 'TPFIL TINE PROFILE') 

IFCILP. EO.l) URITE(3,8) C IFILEC I ), 1-1,29 ) 

8 F0RHATCT2,29A1 ) 

IFCILP. EG. 1) URITE(3,9) 

9 FORMAT CT7, 'M0',T11, 'DAY',T15, 'HR',T19, 'HIN',T23, 'SEC' ) 

C tttttttttt SECTION B XXXXXXXXXXXXXXXXXXXXXXXXtXXX 

CXXXXXXXXXX CHECK TAPE DATE AND SITE SUBTYPE 

c XXXXXXXXXX THE NEXT 8 STATEMENTS FINDS INPUT DATE ON TAPE 

10 CALL TGETCLUN, ISTAT,DATAC 1 ),6510 ) 

CALL TUAITCLUN, ISTAT,NUORDS, IRC) 

IFC IRC. EO.l) GO TO 900 
ICHKCD-DATAC 12819) 

ICHK(2)-DATAC 12820) 

rail t Tf-nr 1 1 

IFCICHKCl ) . GE. IDATEC5 ) . AND . ICHKC2 ) .GE. IDATE C6 ) ) GO TO 20 
GO TO 10 

C XXXXXXXXXX NOW CHECK DATA BLOCK FOR INPUT DATE AND 

C XXXXXXXXXX CALCULATE THE STARTING RECORD POSITION 

20 NRC-1 

J*1 

23 DO 25 1-1,210 

DATAR1 ( I ) -DATA ( J ) 

25 J-J+l 

CALL EBCASCCDATARK1 ),18) 

CALL CFP360CDAT1 < 1 ),48,0 ) 

IFC DATAR1 C 9 ) . EG . IDATE C 5 ) . AND . DATAR1 ( 10 ) . EG . IDATE ( 6 ) ) GO TO 40 
IFC J.GE. 13020) NRC-1 
IFCJ.GE. 13020) GO TO 10 
NRC-NRC+l 


GO TO 23 

C tttttttttt NOU DIRECT CONTROL TO SITE SUBTYPE DATA RETRIEVAL 

40 NSUB*DATAR1 ( 1? >-48 

NSUB1 • DATAR1 C 17 ) 

go to ( 50,50,60,60,60,60 ) nsub 

cmmtm section c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxt t 

c xxxxxxxxxx data retreival for sdas subtype i xmmmmtmmmu 

C ttXXXXXXtt CALCULATE STARTING POSITION FOR RECORD(NRC) 

50 NSTRT-1 

NSTRT-NSTRT+ C ( NRC*210 )-210 ) 

NSTOP-NSTRT+209 

Cttttttxttx INSERT ACTUAL RECORD IN OPERATING RECORD ( DAT AR1 ) 

J«1 

53 DO 55 I«NSTRT,NSTOP 
DATARK J)*DATA(I ) 

55 J-J41 

C XXXXXXXXXX INSERT DATAR1 IN DATAR FOR USE UITH TIRE CONVERSION ( T INCUT ) 

DO 56 1*1,18 

56 DATAR( I )=DATAR1 ( I ) 

C XXXXXXXXXX CONUERT TIRE RECORD AND DATA RECORD 

CALL EBCASC( DATAR ( 1 ), 18 ) 

CALL EBCASCCDATARl ( 1 ), 18 ) 

C XXXXXXXXXX CONUERT EBCDIC UORDS TO PDP WORDS 

CALL CFP360 ( DAT1 ( 1 ) 48 0) 

Cttttttxttx CHECK DATE^ON^TAPEtDATARlO AND 10)) AGAINST INPUT DATE ( IDATEC5 A 
ND 6)) 

IF(DATAR1(9 ) «NE. IDATE(5 ) .0R.DATAS1 ( 10 ) .NE. IDATEt 6 ) ) GO TO 500 
C tttttXtttt CONUERT TINE FROPI ASCII CHARACTERS TO DECIMAL EQUIVALENT 

CALL TinCUT(TIHE) 

Ctttttttttt NEXT 3 STATEMENTS CHECK FOR ERRORS IN TIRE CONVERSION 
IF(NORECS.EQ.O) GO TO 57 
CHK-TIHE-TSAV 
IF(CHK) 303,330,57 

cmtmui output converted record to systems disk for storage 

57 URITE(2'NRDAT,ERR*?O0) TIME, (DATKI), 1-1,48) 

Cttxtttrtxt LINE PRINTER OUTPUT CHECK 

D*FL0AT<NSTRT ) 

IFdLP.EQ. 1 ) URITEC3,59) (DATARKI ), 1*7,16), NRC,D 


CXXXXXXXXXX SITt UATA Tint KKUHLt UUIKU1 IU SUKttN 
YS-YS-15. 

CALL TPLOT(XM,YS,O,0> 

URITEC 1,69 ) CDRTARKI ),I"7,16) 

CXXXXXXXXXX NEXT 5 STATEMENTS CHECK SCREEN LOCATION FOR PROPER OUTPUT 
IF (YS.LT.45. ) XH-XM+150. 

IF( YS.LT.45. ) JRK-JRK+1 
IF(XM.GT.90O. ) XM-50. 

IF (YS.LT.45. ) YS-YSAU 

IF( YS.LT.45.. AND. JRK.GT. 6) YS-700. 

CXXXXXXXXXX INCREMENT THE NUMBER OF STORED RECORDS ( NORECS ) 

NORECS -NORECS* 1 

59 F0RMAT(T2, 'TIME* '5(2A1,2X),3X, 'RECORD NO. ■ 7 ,3X, 14, 3X, 'START POSITI 
XON *',3X,F8.0) 

CXXXXXXXXXX CHECK FOR BLOCK DATA USED UP,IF(NRC.GT.62) GET ANOTHER BLOCK 
IF(NRC.GE .62 ) GO TO 99 

CXXXXXXXXXX INCREMENT THE RECORD NUMBER (NRC) 

NRC-NRC+1 

TSAU-TIME 

CXXXXXXXXXX 270 RECORDS IN 1 DAY, GREATER THAN THAT MEANS STOP 
IF ( NORECS . GE . 270 ) GO TO 500 

CXXXXXXXXXX RETURN AND CALCULATE DATA RECORD STARTING POSITION (NRC) 

GO TO 50 

CXXXXXXXXXX NEXT 5 STATEMENTS RETRIEVE ANOTHER DATA BLOCK 
CXXXXXXXXXX AND RESETS RECORD NUMBER TO 1 
99 CALL TGET(LUN, ISTAT,DATA(1 ),6510) 

CALL TUAIT ( LUN, ISTAT, NU0RD5, IRC ) 

IF ( IRC. EG. 1 ) GO TO 900 

NRC-1 

GO TO 50 

CXXXXXXXXXX SECTION D xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
CXXXXXXXXX DATA RETREIUAL FOR SDAS SUBTYPES 3, 4, 5, 6 XXXXXXXXXXXXXXXXXX 
CXXXXXXXXXX CALCULATE STARTING POSITION FOR RECORD (NRC ) 

60 NSTRT-1 

61 NSTRT-NSTRT*((NRCX210>-210> 

NSTOP-NSTRT+209 

CXXXXXXXXXX INSERT ACTUAL RECORD IN OPERATING RECORD 
J-l 


63 


DO 65 I *NSTRT , NSTOP 
DATARl(J)-DATAd) 

IFd.LT. 12811) DATARS(J)-DATA(I+21d) 


65 J-J+l 

Ctttttttttt CONVERT OPERATING RECORD FROM EBCDIC CHARACTERS TO ASCII CHARACTER 
S 


CALL EBCASC(DATAR1 (1 )*18) 

IFtNSTRT.LT. 12811) CALL EBCASC(DATAR2U ),18) 

C tttttttttt CONVERT EBCDIC DATA UORDS TO PDF WORDS 

CALL CFP360CDAT1 ( l ),48,0) 

IF(NSTRT.LT. 12811) CALL CFP360 CDAT2( 1 ),48,0) 
IF(NSTRT.GE. 12811) GO TO S3 
Ctttttttttt CHECK FOR CONTIGUOUS FiO'333 
IF(DATAR1(17).EQ.DATAR2^ 1?>) S3 70 73 
CmtStmt CHECK DATE )N TO?* AGAINST INPUT DOTE 
IF ( DATAR2 ( 9 ) . NE . I DATE C 5 ) ♦ SR. DAT AR2 < 10 ) . UE . ID ATI 
IF(DATAR1(17).£G.S3) G3 73 S3 


( 3 )) 


GO TO 500 


Ctttttttttt INSERT D ATARI IN DA 70.1 FOR USE UITIi TIDE GQTivER5I0N(TIKCUT ) 

DO 66 1*1,18 

66 DATARC I )*DATAR1 ( I ) 

Ctttttttttt CONVERT TIRE FROM ASCII CHARACTERS TO DECIRAL EG’J I VALENT 

CALL TlflCVT(TinEl) 

Ctttttttttt INSERT DATAR2 IM DATA? ~CR USE UITH TIRE CONVERSION(TIRCVT) 

DO 67 1*1,18 

67 DATAR(I )*DATAR2(I ) 

cmmtftt CONVERT TIRE FROM ASCII CHARACTERS TO DECIMAL EQUIVALENT 
CALL TIHCVT(TIflE) 

NK-NCH-48 

Ctttttttttt CHECK FOR TIRE SEQUENCE 

IFCTIHE1.NE.TIHE ) GO TO 403 

Cttmmtt NEXT 3 3TATERENT3 CHECK FOR ERRORS IN TIKE ONVERSIOH 

IF(NORECS.EO.O) GO TO 63 

CHK-TIHE-TSAV 

IF(CKK) 800,300,63 

Ctttttttttt OUTPUT CONVERTED RECORD TO S73TERS DISK FOR STORAGE 

68 URITE ( 2 ' NRDAT , ERR *700 ) TIRE, (DATld ), 1*1,43), CDA72d ), 1*1, NK) 

Ctttttttttt LIME PRINTER OUTPUT CRICX 

D-FLOAT(NSTRT) 


IF(ILP.EQ.l) URITE(3,59) (DATAR1CI >, 1-7,16), NRC,D 
C tttttttttt SITE DATA TIHE PROFILE OUTPUT TO SCREEN 

VS-YS-1S. 

CALL TPLOT(XH,YS,e,0> 

URITE( 1,69) (DATAR1 (I ), 1-7, 16 ) 

C tttttttttt NEXT S STATEMENTS CHECK SCREEN LOCATION FOR PROPER OUTPUT 

69 F0RMAT(1H+,1( (2A1, '/' >,2A1,2X,2(2A1, ' i ' ),2A1,2X) > 

IF(YS.LT.45. ) XM-XH+150. 

IF(XH.GT .900* ) XM-50. 

IF(VS.LT.45. ) JRK-JRK+1 
IF(YS.LT.45. ) YS-YSAV 
IF(YS.LT.45. .AND. JRK.GT.6) YS-730. 

C tttttttttt INCREMENT THE NUMBER OF OUTPUT RECORDS(NORECS ) 

NORECS-NORECS+1 
TSAV-TIME 


C tttttttttt CHECK FOR BLOCK DATA USED UP,IF(NRC.GE.61 ) GET ANOTHER BLOCK 

IF(NRC.GE.61) GO TO 73 

C xtxtxtttxt INCREMENT THE RECORD NUMBER (NRC) 

NRC-NRC+2 

C tttttttttt 270 RECORDS IN 1 DAY, GREATER THAN THAT MEANS STOP 

IF ( NORECS .GE .270 ) GO TO 500 

ctmmm return and calculate data record starting position (nro 

GO TO 60 

c tttttttttt NEXT 5 STATEMENTS RETRIEVE ANOTHER DATA BLOCK 

c tttttttttt AND RESETS RECORD NUMBER TO 1 

73 CALL TGET(LUN,ISTAT,DATA(1),6510) 

CALL TUAIT( LUN, ISTAT, NUORDS, IRC ) 

IF ( IRC .EG. 1 ) GO TO 900 
NRC-1 
GO TO 60 

C tttttttttt CONTROL TRANSFER TO HANDLE CONTIGUOUS RECORDS 

75 Nl-NSTRT+16 

76 Nl-Nl+216 
CALL EBCASC(DATA(N1 ),18) 

IF(DATACNl).EG.NSUBl) GO TO 76 
NS- C CN1-16 )-NSTRT )/21t 
NRC-NRC+(NS-1) 

YS-YS-15* 


c c 


g5 

<o 

5 Q 

^ t f'i’f 


CALL TPLOT ( XS , YS, 0, 0 ) 

WRITE (1,77) NS 

77 FORHATC 1H+, 14, SX, 'CONTIGUOUS RECORDS UERE ENCOUNTERED' ) 

IF(NRC*LE»62 ) GO TO 60 

CALL TGET(LUN,ISTAT,DATA(1),6510) 

CALL TUAIT(LUN, ISTAT,NUORDS, IRC ) 

IF(IRC.EQ.l) GO TO 900 

NRC-1 

GO TO 60 

C tttttttttt CONTROL TRANSFER TO HANDLE A SUBTYPE SHIFT IN BLOCK 

80 NSTRT*211 
GO TO 61 

C ********** CONTROL TRANSFER TO HANDLE SUBTYPE 3 AND 5 AS LAST RECORD IN BLOCK 

90 CALL EBCASC(DATAR1(1 ), 18) 

CALL CFP360(DAT1 ( 1 ), 43, 0 ) 

DO 95 1-1,18 
95 DATARC I )*DATAR1 ( I ) 

CALL TINCUT(TIflEl) 

CALL TGET(LUN,ISTAT,DATA(1 ),S310) 

CALL TUAIT ( LUN, ISTAT, NUORDS, I RC ) 

IF(IRC.EQ.l) GO TO 900 
DO 100 1-1,210 
100 DATAR2< I )-DATA(I ) 

CALL EBCASC(DATAR2( 1 ), 13 ) 

CALL CFP360(DAT2( 1 ),43,3) 

DO 105 1-1,18 
105 DATAR(I )*DATAR2( I ) 

CALL TINCUT(TinE) 

IF(TIHEl.NE.TIME) GO TO 400 

URITE ( 2 ' NRDAT, ERR-700 ) TINE, (DAT1 (I ), 1-1,43 ), (DAT2C I ), I»1,NK ) 

YS-YS-15. 

CALL TPLOT(XN,YS.0,0) 

URITE ( 1,69 ) (DATARl(I), 1-7,16) 

NORECS-NORECS+1 

NRC-1 

NSTRT-211 

IF ( NORECS • GE • 270 ) GO TO 500 


60 TO 61 

Ctttttttttt ERROR MESSAGES AND OUTPUT INFO BELOUttSSSUS*ttX»S«ttStt» 
400 YS-YIC-15. 

CALL TPLOT ( XS, VS, 0, 0 ) 

URITE( 1,405) (DATARKI ),I-?,16) 

40S FORMATdH*, 'THIS DATA IS OUT OF TINE SEQUENCE AT', EX, 

X3C2A1, '✓' ),3X,3(2A1, ' * ' ) ) 

GO TO 1000 
S00 YS-YK-15. 

CALL TPLOT CXS, VS, 0,0) 

WRITE (1,505) NORECS 

505 FORMAT (1H+, 'THE DATA HAS BEEN STORED WITH ', 14, 2X, 'RECORDS' ) 

GO TO 1000 
600 YS-YIC-15. 

CALL TPLOT(XS,YS,0,0) 

WRITE (1,605) 

60S F0RMAT(1H+, 'THERE IS AN ERROR IN FILE MANIPULATION') 

GO TO 990 
700 YS-YK-15. 

CALL TPLOT (XS,YS,0,0) 

WRITE (1,705) 

705 FORMAT ( 1H+ , ' THERE IS AN ERROR IN WRITING THE FILE') 

GO TO 1000 
800 YS-YK-15. 

CALL TPLOT(XS,YS,O,0) 

WRITE(1 805) 

805 FORHATUH+, 'THERE IS AN ERROR IN THE TIME CONVERSION' ) 

GO TO 1000 
900 YS-YIC-15. 

CALL TPLOT (XS,YS, 0,0) 

WRITE (1,905) 

HFILE-HFILE+1 

905 FORMAT CIH*. 'AN END OF FILE WAS ENCOUNTERED DURING EXECUTION') 

1000 VS- VS- 15. 

CALL TPL0T(X$,VS,0,0) 

CALL TPAUSE 
CALL ERASE 
WRITEC 1,2000) 


99e 

leei 
2000 
2050 
9999 
PI P> 


READ (1,2050) IREP 

cEosfjSiM^gSiSg^ 3,1 3ISP0SE- ^ ' ?RI * T '‘ •' ESR -" 99> 

IF(IRcP.E0.1 ) GO TO 2 
CALL RUND ( ISTAT ( 1 ) ) 

URITE(1,1S01) 

F0RHAT(/V1X, 'THE TAPE IS BEING REWOUND' ) 

V0U UISH TO RETRIEVE MORE SATA ? 


FORflATdl ) 
STOP 


d*VES,2«N0)' ) 


END 


4 *. 


SUBROUTINE TAPFUD(NFILE,FIFILE,ISTAT,IREP) 
DIMENSION ISTAT(l) 

LUN-7 

IF C IREP . EG . 2 . AND . NFI LE . EQ ♦ 1 ) RETURN 
IF(NFILE.EQ.l) IDIR-3 
IF(NFILE*EQ.l ) GO TO 60 
IFCNFILE-nFILE) 20,20,50 
20 NF*MFILE-NFILE+1 
IDIR-13 
GO TO 60 

50 NF-NFILE-PIFILE 
IDIR-10 

60 CALI SFUN ( LUN, I$TAT,NF, IRC, IDIR ) 

CALL TUAIT( LUN, ISTAT,NUORDS, IRC ) 

RETURN 

END 


SUBROUTINE TINCUTCTINE) 
INTEGERS^ ITIBE(3),ITINE1(6) 
LOGICALS1 DATAR(210) 

CORHON/T I R/D AT AR 
IC-11 

5 H-IC-IO 

DO 10 J-48,S? 
IF(DATAR(K).EQ.J)ICHK"J 
10 CONTINUE 

ITIHE1 (H )*ICHIC-48 
IFtH.GT .6) GO TO IS 
K»K+1 
GO TO S 
IS J-l 

DO 20 1*1,5, 2 
ITIBEKI )*ITINE1(I )*10 
ITIHE( J )*ITINE1 C I HITinil (1+1 ) 
J»J+1 

20 CONTINUE 

Tl*(FL0ATCITINE(3))/63. ) 

T2« (FLOAT ( ITINE(2 ) )+Tl )/60. 
TINE* (FLOATC ITIBEC 1 ) HT2 ) 
RETURN 
PIP> END 


TAPFIL VARIABLES 


CHK- 

D- 

IDIR- 

ILP- 

IRC- 

IREP- 

IRK- 

LUN- 

MFILE- 

Nl- 

NCH- 

NF- 

NFILE- 

NK- 

NORECS 

NRC- 

NRDAT- 


NS- 

NSTOP- 

NSTRT- 

NSUB- 

NSUBl- 


used for the difference between two consecutive time points 

real value of starting position in data block for conversion and storage 

tape function to be performed (3« rewind, 10* forward space flies, 

13* backspace files) 

Integer variable for line printer output decision (1- yes, 2" no) 

denotes error status (0« correct read, 1« end of file, 2“ tape error) 

determines whether want more data (l»yes, 2" no) 

denotes the screen time output column number (1-6) 

logical unit number of tape drive 

current file number 

subscripted variable designating subtype value; used to check for 
contiguous records 

the number of channels at a site 

number of files to space on the IBM 360 tape 

denotes the tape file number in which the desired date Is located 

the number of channels of data contained In subtype records 4 and 6 

the number of output records on the system disk 

the operating record number in the data block from the tape 

associate variable for opening direct access file (JFILE) and for 
outputting to the system disk 

the number of contiguous records 

stopping position in data block, used with NSTRT (NSTOP*NSTRT + 209) 

starting position in data block for conversion and storage; 
calculated from NRC 

denotes the number of the subtype 

the ASCII code equivalent of the subtype 


50 


NT- 


the number of channels + 1 for time 


NWORDS- the number of words not read from the block on the tape 

Tl- converts excess seconds to minutes in TIMCVT 

T2- converts excess minutes to hours in TIMCVT 

TIME- decimal equivalent time value converted from ASCII by TIMCVT 

TIME1- decimal equivalent time value converted from ASCII characters stored In 

subtypes 4 and 6 

TSAV- used to save previous time point for use with CHK to check time conversions 
XM- T4014 scope column abscissa value for time output (50-800) 

XS- T4014 scope abscissa location used for output 

YK- T4014 scope ordinate value used at end of execution for output (YK*40) 

YS- T4014 scope ordinate location used for output 

YSAV- T4014 scope ordinate value used to save last value of YS 


DAT1(48)- 

DAT2(48)- 

DATA( 13020)- 

DATAR(210)- 

DATAR1(210)- 

DATAR2(210)- 

ICHK(2)- 

IDATE(6)- 

IFILE(30)- 

ISTAT(4)- 

ITIME(3)- 

ITIME1(6)- 


TAPFIL ARRAYS 

real storage of subtype 1,3, and 5 records equlvalenced to 
DATARl (19-210) (channels 1-48) 

real storage of subtype 4 and 6 records equlvalenced to 
DATAR2 (19-210) (channels 49-96) 

LOGICAL *1 array to store IBM 360 words from tape 

LOGICAL *1 array for time conversion in subroutine TIMCVT 

LOGICAL *1 array to store individual subtype 1, 3, and 5 records 
(channels 1-48) 

LOGICAL *1 array to store individual subtype 4 and 6 records 
(channels 49-96) 

LOGICAL *1 array used to locate user date on tape 

LOGICAL *1 array of desired date to be retrieved (Input from scope) 

LOGICAL *1 array used to store data file name for storage 

status code for tape handling 

stores decimal equivalent of hours, minutes, and seconds 

stores decimal equivalent of each character of desired time 

51 


APPENDIX D 

CHPLOT SOURCE LISTING 


PRECEDING PAGE BLANK NOT FILMED 
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SUBROUTINE AU70SC ( XHAX, XMIN, STX, RLX, TICK ) 
INTEGERS4 IA,ISTX.IB,IEX 

XHAX (INPUT ARGUMENT ) MAX. UALUE 

XHIN (INPUT ARGUMENT) MIN. UALUE 

STX (RETURN ARGUMENT) STARTING X UALUE ON AXIS 

RLX (RETURN ARGUMENT) LENGTH OF AXIS 

TICX (RETURN ARGUMENT) LENGTH BETWEEN TIC MARKS 

RLX-XMAX-XHIN 

QaRLX 

IF(RLX.EQ.O.) GO TO 1 
IA-ALOGIO(RLX) 

IF(IA)1,1,3 

CONTINUE 

ISTX-XHIN/(10.«IA) 

IB-ISTX 

IF(XMIN.LT.0. )IB«ISTX-1 
ISTX-IB*(10UIA) 

IEX a XMAX/( lO.StlA ) 

IB-IEX+1 

IF(XMAX.LT.e.)IB-IEX 

IEX»IB*(10«IA) 

RLX-IEX-ISTX 

STX-ISTX 

TICX» ( IEX-ISTX )/10 
IF(G.LE.10. )TICX-1. 

GO TO 4 
CONTINUE 

IF(RLX.GE.l. )GO TO 3 
IA-IA-1 

ISTX»XHIN/(10.t*IA) 

IB-ISTX 

IF(XMIN.LT.O. )IB*ISTX-l 
STX-IBt(lO.KIA) 

IEX-XHAX/(10.»IA) 

IB-IEX+1 

IF(XHAX.LT.». )IB-IEX 
EX-IBt(lB.UIA) 


RLX-EX-STX 

TICX*(i0.«IA> 

8 IF(TICX.GE.<RLX/10. 

TICX-TICX+C10.MIA) 
GG TO 8 
7 CONTINUE 

IF C (RLX/TICX ) . GT . 4 . 
TICX-TICX/3. 

GO TO 7 
4 RETURN 

END 

PIP> 


>>GO T07 
)GO TO 4 


CHPLOT: 


DIMENSION SX 1 ( 6 ) , SX2 ( 6 ) , SX3 ( 6 ) , SX4 ( 6 ) , SXS ( 6 ) 

DIMENSION SAW 1 ( 300 ) , b AU2 ( 300 ) , SAU3 ( 300 ) » SAU4 ( 300 ) , SAVE ( 300 ) 
DIMENSION PID(4),XTITL(4),TIHE<300) 

C0MH0N/SX/SX1 , SX2 , SX3, SX4, SX5 

COMNON/PIDi/PID 

COMHON/STOPT/ISTOPT , ISTOPF 

COMMON/POINTS/NPT 

COHHON/HOUM/NPLT, IFILE,NSCH 

COMHON/LUN/IS, IP, IC, ISD 

COHHON/UTL 1 / IFLAGX, XH AX , XU I N 

COMHON/UTL2/ IFLAGY, YHAX , YM I N 

C0HH0N/NUMBR/SAU1 , SAU2 , SAU3, SAU4 , SAU5 , TIME 

DATA XTITL/'CLOC'.'IC TIS'ME S' V 

CALL BEGIN (9600) 

IS-1 

IST0PT*0 
ISTOPF-O 
1 CALL DISPLY 
IFLAGX-0 
IFLAGY-0 
CALL RETREU 

5 CALL PLOTS < PI D (1 ) , *TI TL ( 1 ) , S AV 1 ( 1 ) , SAU2 ( 1 ) , S AU3 ( 1 ) , SAU4 C 1 > , SAU5 ( 1 ) 
*.TINE< 1 ) , NPT , NPT , NPT, NPT , NPT, - 1 ) 

CALL TPAUSE 
CALL UATNXT(IR) 

ISTOPT-1 
ISTOPF- 1 

IFdR.EQ.l ) GO TO 1 

ISTOPT-O 

ISTOPF-O 

IF(IR.EQ.3) GO TO 1 
ISTOPF -l 
ISTOPT-O 

IF(IR.E0.4> GO TO 1 
IFCIR.E0.5) GO TO 10 
CALL PSCALE 

GO TO 5 
10 STOP 
END 


cm 

cm 

cm 

cm 

cttt 


SUBROUTINE CHSIZE(N) 

THIS ROUTINE CHANGES CHARACTER SIZE ON A TEKTRONIX 4014-1 DISPLAY 
SCOPE ACCORDING TO N AS FOLLOWS: N-l - LARGEST 

N-2 - NEXT SMALLER 
N-3 - NEXT SMALLER 
N»4 - SMALLEST 

NP-N+55 
CALL CHOUT ( 2? ) 

CALL CHOUT (NP) 

RETURN 

END 


PIP> 


SUBROUTINE DATIN(XS,VS,DATU) 

CXtt THIS ROUTINE ALIOUS THE USER TO RETRIVE FLOATING POINT DATA FROH 

CXXX THE SCOPE IN G13.6 FORHAT, ONE UORD PER CALL. XS IS THE ABSOLUTE 

CXXX SCREEN COORDINATE IN TEKPOINTS, LIKEWISE VS. DATA FROH THE SCOPE 
CXXX IS RETURNED IN DATU. 

COHHON/LUN/IS, IP, IC, ISD 
CALL SCALEU. 0,1. 0,0. 0,0.0) 

CALL ALPHA 

CALL TPLOT(XS, YS,0,0 ) 

READ (IS, 100) DATU 
100 F0RNAT(G13.6) 

RETURN 

END 


PIP> 


i 


i 

! 


SUBROUTINE DISPLY 

DINENSION SX1 ( 6 ) , SX2 ( 6 ) , SX3 ( 6 ) , SX4 ( 6 ) , SX5 ( 6 ) 

DIMENSION PID<4 ) 

LOGICALJ1 IFILE(30), IFILET(SO ) 

COHNON/CHNO/ NCH1 , NCH2, NCH3, NCH4, NCH5 
COHHON/PID1/PID 
COMMON/STOPT/ISTOPT, ISTOPF 
CONMON/SX/SX1 ,SX2,SX3, SX4, SX5 
COMMON/HOUM/NPLT, IFILE.NSCH 
COMHON/CHECIC/ICHECIC.ANPLT 
CONHON/FILEB/IFILET 
DATA IFILE(30)/0/ 

DATA IFIlET(50>/0/ 

CALL ERASE 

CALL SCALE<1.0,1. 0,0. 0,0,0) 

CALL ALPHA 
CALL CHSIZE(2 ) 

CALL TPLOT ( 25 . , 650 .,0,0) 

URITEU.S) 

5 F0RMAT(1H+, 'EL 55 CHANNEL PLOTTING PROGRAM' ) 

CALL CHSIZEO ) 

XS-20. 

YS-650. 

IFdSTOPF .EQ.l ) GO TO 16 
YS-YS-50. 

CALL TPLOT(XS, YS,0,0 ) 

WRITE (1.10) 

10 F0RMAT(1H+, "ENTER PLOT FILE NAME ' ) 

YS-YS-25. 

CALL TPLOT (XS.VS.0.0) 

READC1.15) (IFILECI). 1-1,29) 

15 FORMAT (29A1) 

16 YS-YS-25. 

CALL TPLOT(XS, YS.0,0) 

WRITE (1,20) 

20 F0RHAT(1H+, 'ENTER NO. OF PLOTS ON THIS GRAPH (MAX-5)' ) 
YS-YS-25. 



CALL TPLOT(XS,YS,0,O) 

READ(1,25) HPLT 
25 F0RHATC12) 

ANPLT-FLOAT(NPLT) 

IF(ISTOPT.EQ.l) GO TO 29 
yg.yS.25 , 

CALL TPLOT(XS,YS,O,0) 

URI TEC 1,26) 

26 FORMAT ( 1H+, 'ENTER PLOT CAPTION (49 CHARACTERS ) ' ) 
YS-YS-25. 

CALL TPLOT (XS,YS,0,0) 

READ (1,2?) CIFILETCI), 1-1,49) 

27 FORMAT (49A1) 

YS-YS-25. 

CALL TPLOT(XS,YS,0,O) 

URITE(1,43) 

43 FORMAT (1H+, 'ENTER CAPTION FOR LEFT SIDE OF PLOT * ) 
YS-YS-25. 

CALL TPLOT (XS,YS, 0,0) 

READ (1,44) (PID(I), 1-1,4) 

44 FORMAT (4A4) 

29 YS-YS-25. 

CALL TPLOT(XS,YS,O,0) 

WRITE (1,30) 

30 FORMAT ( 1H+, 'ENTER CHANNEL NO. AND UNITS' ) 

IC-1 

RS-XS+40. 

35 YS-YS-25. 

GO TO (40,45,50,55,60) K 

40 CALL TPLOT (XS, VS, 0,0) 

READ(1,41> NCH1 

41 F0RHATU2) 

CALL TPLOT (RS, VS, 0 , 0 ) 

42 FORMAT (6A4) 

READ (1,42) (SX1CI), 1-1,6) 

IF(K.EQ.NPLT) GO TO 70 

K-K41 

GO TO 35 






45 CALL TPLOTCXS, YS,0,0) 

READ (1,41) NCH2 
CALL TPLOTCRS, YS, 0,0 ) 

READ (1,42) (SX2CI ), 1*1,6 ) 

IFCK.EQ.NPLT)GO TO 70 

K-K+l 

GO TO 35 

50 CALL TPLOTCXS, YS, 0,0 ) 

READC 1,41 ) NCH3 
CALL TPLOT(RS,YS,0,0) 

READC 1,42 ) (SX3(I), 1-1,6) 

IFCK.EQ.NPLT ) GO TO 70 

K-K+l 

GO TO 35 

55 CALL TPLOT(XS,YS,0,0) 

READC 1,41 ) NCH4 
CALL TPLOTCRS, YS,0,0 ) 

READ (1,42) (SX4C I ),I*1,6) 

IFCK.EQ.NPLT) GO TO 70 

K-K+l 

GO TO 35 

60 CALL TPLOT(XS,YS,0,0 ) 

READC 1,41) NCH5 
CALL TPLOTCRS, YS, 0,0) 

READC 1,42 ) CSX5CI), 1-1,6) 

70 YS-YS-25. 

CALL TPLOTCXS, YS,0,0 ) 

IFC ISTOPF.EQ. 1 ) GO TO 75 
URITEC1,71 ) 

71 FORNAT C 1H+, 'ENTER NO. OF CHANNELS AT SITE') 
YS-YS-25. 

CALL TPLOTCXS,YS,0,0) 

READC 1,41) NSCH 
VS-vs-30. 

CALL TPLOTCXS, VS,0,0) 

V~ URITEC 1,72 ) 

72 FORHATC1H4, 'DATA IS BEING RETRIEVED' ) 

RETURN 

END 


1 

: 

■| 

i 



SUBROUTINE PLOTS(PID,UNITS,Yl, Y2,Y3,V4,YS,X,NPTl,NPT2,NPT3,NPT4, 

XNPT5, IG ) 

cm PL0T6 BUILDS 5 PLOTS OF DATA ON A STANDARA GRAPH FOR USE IN THE 

Cttt QUICK-LOOK GRAPHICS PROGRAM. PID IDENTIFIES THE NAME OF THE DATA BASE, 

cm UNITS IS THE ARRAY CONTAINING THE UNITS OF MEASUREMENT OF THE 

CSX* DEPENDENT PARAMETER. XTITL IS THE ARRAY CONTAINING THE INDEPEDENT 

CXXX PARAMETER UNITS. X AND Y ARE THE ARRAY ELEMENTS FROM THE FIRST ARRAY 
CXXX OF DATA READ IN, CORRESPONDING UITH SYMBOL XKENO. X2 AND Y2 ARE THE ARRAY 
CXXX ELEMENTS FROM THE SECOND ARRAY OF DATA, CORRESPONDING UITH SYMBOL XKVNO. 
CXXX NPID IS THE NUMBER OF CHARACTERS TO BE PRINTED FOR PID. NUN IS THE NUMBER 

CXXX OF CHARACTERS TO BE PRINTED FOR UNITS. NXT IS THE NUMBER OF CHARACTERS 

CXXX TO BE PRINTED FOR XTITL. NPXY IS THE NUMBER OF POINTS IN THE FIRST 

CXXX DATA ARRAY, NPXY2 IS THE NUMBER OF POINTS IN THE SECOND DATA ARRAY. 

CXXX IG IS THE GRID LINE VARIABLE, (IG. EG. -1 )FOR GRID LINES. (IG.EQ.l) 

CXXX FOR NO GRID LINES. 

DIMENSION PID(4),UNITS(4),Y1(1 ),X(1 ),Y2(1 ),Y3(1 ) 

DIMENSION SX1(6),SX2(6),SX3(6),SX4(6),SX5(6) 

DIMENSION FILET (12) 

DIMENSION Y4(1),Y5(1> 

COHMON/F I LEB/ IF I LET 
COHHON/SX/SX 1 , SX2, SX3, SX4 , SX5 
COMMON/CHECK/ ICHECK , ANPLT 
COMMON/LUN/IS, IP, IC, ISO 
LOGICALX1 IFILET(50) 

COMMON/UECT/IUECT 
EQUIVALENCE (FILETt 1 ),IFILET( 1 ) ) 

DATA IFILET(S0)/O/ 

IQ IVECT-1 

CALL SSPLOT(Yl(l ),V2(1),Y3(1 ), Y4(l ),Y5(1 ),X(1 ),850., 

X623. , 140., 1 18. ,NPTl,NPT2,NPT3,NPT4,NPT5, IG, IG ) 

CALL CHSIZE(3) 

CALL VLABEL(10.,?00.,PID(1),16) 

CALL XLABEL(420.,S8.,UNITS(1),16) 

CALL SCALE(1.0,1.0,0.0,0.0) 

CALL TPLOT(10.,10.,0,0) 

CALL ALPHA 


UK 1 TUI, by) UHLtT(i),l-l,4y) 
59 FORHAT ( 1H*# 49AI ) 

NPLT-INT(ANPLT) 

DO 100 J-1,NPLT 
GO TO (101,102, 103, 104,105) J 
101 CALL TPL0TC634. ,70. ,0,0) 

CALL VECTOR 

CALL TPLOT(710.,70.,1,0) 

CALL POINT 

aai t At PUA 

CALL XLABEL(753.,70.,SX1(1 ),24 ) 
GO TO 100 

102 CALL TPLOT ( 634 .,50 .,0,0) 

CALL VECTOR 
CALL CHOUT ( 27 ) 

CALL CHOUT (97) 

CALL TPLOT(710.,50.,1,0) 

CALL POINT 
CALL ALPHA 

CALL XLABEL ( 753 . , 50 • , SX2 ( 1 ) , 24 ) 
GO TO 100 

103 CALL TPLOT (634* ,30. ,0,0 ) 

CALL VECTOR 
CALL CHOUT (27) 

CALL CHOUT (99) 

CALL TPLOT(?10.,30.,1,0) 

CALL POINT 
CALL ALPHA 

CALL XLABEL(?S3.,30.,SX3(1),24> 
GO TO 100 

104 CALL TPLOTC113.,70.,0,0) 

CALL VECTOR 
CALL CHOUT (27) 

CALL CHOUT (100) 

CALL TPLOTC194.,70.,1,O> 

CALL POINT 
CALL ALPHA 

CALL XLABELC237.,?0.,SX4C1),24) 



W TV 1V« 

1*5 CALL TPLOT(il8»*SO* # O # 0) 

CALL VECTOR 
CALL CH0UTC27) 

CALL CH0UTC98) 

CALL TPLOTC194. # SO..1,0> 

CALL TPLOTC194. # 5O..0,O> 

CALL POINT 

CAf f ALPHA 

CALL XLABEL(237.,S0.,SXS<1),24) 
189 CONTINUE 

CALL TPLOT<237.,50.,0.0> 

CALL ALPHA 
CALL HONE 
RETURN 
PIP> END 


SUBROUTINE PSCALE 

COHNON/UTL 1 / IFL AGX , XN AX , XN IN 

COMMON /UT L2 / IFLAGY,YHAX,YMIN 

CALL ERASE 

IFLAGX-1 

IFLAGY-1 

XS-2S. 

VS -600. 

CALL TPLGT<XS,YS,0,0) 

WRITE! 1 20) 

20 FORMAT! IH+, 'ENTER NAXINUN AND NININUN UALUES FOR X AND Y') 
YS-YS-25. 

CALL TPLOT ( XS, YS, 0, 0 ) 

WRITE!1,2S) 

25 FORMAT! 1H+, 'XMAX- ' ) 

CALL DATIN(10S.,YS,XHAX) 

CALL ALPHA 

YS-YS-25. 

CALL TPLOTtXS, YS,0,O) 

URITE( 1,26 ) 

26 F0RHAT!1H*,'XMIN-') 

CALL DATIN( 105. , YS,XHIN ) 

CALL ALPHA 

YS-YS-25. 

CALL TPLOT(XS,YS,0,0) 

WRITE (1,87) 

27 FORMAT! 1H+, 'YMAX-' ) 

CALL DATIN( 105. , VS, VNAX ) 

CALL ALPHA 

YS- YS— 25 . 

CALL TPLOTCX$,YS,0,O) 

WRITE (1,28) 

28 FORMAT! 1H4. 'YMIN- ' ) 

CALL DATIN( 105. ,VS,VNIN) 

CALL ALPHA 

RETURN 

END 


SUBROUTINE RETREV 

DIRENSION SAU1C300)#SAU2C300)#SAU3(300)#SAU4C300)»SAUSC300) 
DIRENSION DATAC100)#TIHEC300) 

LOGICALS1 IFILEC30) 

CONNON/NUNBR/SAU1 , SAVE, SAV3# SAU4# SAU5# TINE 
COmWN/HOUn/NPLT, IFILE.NSCH 
CONNON/POINTS'NPT 

COHNON/CHNO/ NCN1#NCH2,NCH3#NCH4,NCH5 
DATA IFILEC30)/0/ 

NT-NSCH+1 

0PENCUNIT*2#NAHE*IFILE#TYPE* 'OLD'# ACCESS* 'DIRECT'# 

{RECORDS I ZE ■ NT # ASSOC I ATE VAR I ABLE ■ NRDAT # ERR* 25 ) 

NRDAT-i 
5 K*1 

10 READ(2'NRDAT#ERR*21 ) (DATACI ), I»1#NT) 

TIRE(IC)*DATA(1 ) 

SAUl(K)*DATA(NCHl4i) 

IFCNPLT.LT. 2) GO TO 20 
SAU2CIC)*DATACNCH2+1 ) 

IFCNPLT.LT. 3) GO TO 20 
SAU3CIC)*DATACNCH3*1 ) 

IFCNPLT.LT. 4) GO TO 20 
SAU4CIC)*DATACNCH4+1) 

IFCNPLT.LT. 5) GO TO 20 
SAUSCK)*DATACNCH5+| ) 

20 IFCIC.LT.270> IC*IC*1 
IFCK.LT.270) GO TO 10 

21 NPT-IC-1 

CLOSE CUNIT-2# ERR-2S > 

2S RETURN 
END 


PIP> 


csss 

csss 

csss 

csss 

csss 

csss 

csss 

csss 


10 

11 

csss 


12 

13 


SUBROUTINE SSPLOTC Y1,Y2, Y3, Y4,Y5,X,TPX, TPY,XO, V0,NPTl,NPT2 
S,NPT3,NPT4,NPT5, IGX, IGY ) 

SSPLOT AUTONATICALLV SCALES THE X AND Y DATA STORED IN ARRAYX AND 
ARRAYY TO FIT INTO THE TEKPOINTS AVAILABLE AS SHOUN IN TPX AND TPY. 
SSPLOT THEN DRAWS THE AXIS AND GRID, PLOTS THE DATA IN ARRAYX 
AND ARRAYY AND LABLES THE NUNERICAL VALUES OF THE TIC NARKS ON THE AXES. 
XO AND YO ARE ORIGIN LOCATIONS FOR THE GRID IN ABSOLUTE SCREEN 
COORDINATES REFERENCED FRON SCREEN (0,0). IGX AND IGY ARE 
INDICATORS FOR GRID LINES ON THE PLOT (.LT.0) OR TIC NARKS 
C.GT.0) ON THE OUTPUT PLOT. 

DIHENS ION VI ( 300 ) , Y2 ( 300 ) , Y3 C 300 ) , Y4 ( 300 ) , Y5 ( 300 ) 

DIHENSION X(300) 

CONNON/LUN/IS, IP, IC, ISO 
COHHON/CHECK/ I CHECK , ANPLT 
COHHON/VECT/ IVECT 
COHHON/UTL1/ IFLAGX,XHAX,XHIN 
C0HN0N/UTL2/ IFLAGY, YNAX, YNIN 
INTEGERS4 N 
CALL ERASE 
L»1 


IF(IFLAGX.GT.0 ) GO TO 11 
XNAX-X(l) 

XHIN-X(l) 

DO 10 I-1,NPT1 

IF(X(I ).GT .XNAX) XNAX-X(I) 

IF(Xd).LT.XNIN) XNIN-X(I) 

CONTINUE 

CALL AUTOSC ( XNAX , XN IN, STX , RLX , TI CX ) 
CHECK PLOTTING EFFICIENCY 
N-RLX/TICX 
DO 12 I»l. N 

IFC (STX-fTICX).GE.XNIN) GO TO 13 

STX«STX*TICX 

RLX-RLX-TICX 

CONTINUE 

CONTINUE 

DO 14 I*1 # N 


ss 


IFC (STX*RLX-TICX).LE.XNAX) 60 TO 16 
RLX-RLX-TICX 
14 CONTINUE 
16 CONTINUE 

INRLX«RLX/TICX 
CRLX-INRLXSTICX 
IFCCRLX.LT.RLX) INRLX-INRLX+l 
RLX*INRLXtTICX+#« 1 
IF ( IFLAGY. GT.0 ) GO TO 21 

C ttttt FIND LIHITING Y VALUES FOR LEFT SIDE StttSSX 

YNAX-VKl) 

VNIN-VlCi) 

NPLT-INTCANPLT) 

DO 17 J-t,HPLT 

GO ?0 (99,101,102,103,104) J 
99 DO 18 I-1,NPT1 

IF(YKI).GT.YNAX) YNAX-YKI) 

IF(YHI).LT.YNIN) YHIN*Y1(I) 

18 CONTINUE 
GO TO 17 

101 DO 19 I»1,NPT2 
IF(Y2(I).GT .YNAX) YNAX-Y2CI) 

IFC Y2 ( I) . LT. YNIN ) YNIN-Y2CI) 

19 CONTINUE 
GO TO 17 

102 DO 20 I-1,N?T3 
IF(Y3(I).GT.YNAX) YNAX-Y3CI) 

IF(Y3(I ).LT.YHIN) YNIN-Y3(I) 

20 CONTINUE 
GO TO 17 

103 DO 200 I-1,NPT4 
IF(Y4(D.CT.VNAX) YFAX-Y4CI) 

IFCY4CD.LT. YNIN) YKIN-Y4CD 

200 CONTINUE 
GO TO 17 

104 DO 201 I-1.NPT5 

IFCY5CI ).GT.VNAX) VNAX-V5CD 
IFCYSCD.LT. YNIN) VNXN-VSCD 


201 CONTINUE 

1? CONTINUE 

21 CALL AUTOSC<YHAX # YNIN,STY,RLY,TICV> 

C*« CHECK PLOTTING EFFICIENCY 

N-RLY/TICY 
DO 22 I-1,N 

IF( (STY+TICYKGE.YHIN ) GO TO 23 

STY-STY+TICY 

RLY-RLY-TICY 

22 CONTINUE 

23 CONTINUE 

DO 25 I*1#N 

IF ( CSTY+RLY-TICV ).LE. YHAX ) GO TO 26 
RLY-RLY-TICY 

25 CONTINUE 

26 CONTINUE 
INRLY-RLV/TICY 
CRLY-INRLYfTICY 
IF(CRLY.LT.RLX) INRLY»INRLY*1 
RLY-INRLY*TICY+0.1 

C ***** CHECK PLOTTING EFFICIENCY ttttttttt 
IF ( CRLVR . LT . RLYR ) INRLYR- INRLYR+1 
RLYR- INRLYR*T ICYR* . 1 
XFACT -TPX/RLX 
YFACT-TPY/RLY 
YOS-YO 
YON-YO 
XOS-XO 

XO-XO-STX*XFACT 
YO-YO-STY*YFACT 
YOR- YOH-STYR*YF ACTR 
IF(YO*GT • (YOS+TPY) ) YO-VOS+TPY 
IFCYOR.GT. ( YOS+TPY) ) YOR-YOS+TPY 
30 CALL SCALE ( XFACT , YFACT, XOS, YOS ) 

CALL AXISC0«0 # 0.0 # RLX # RLY # TICX*TICY # -1,-1 ) 
CALL SCALE C XF ACT # YF ACT, XO • VO ) 

CALL POINT 

IF(IUECT.GT.0) CALL VECTOR 


1-1 

NPLT * I NT < AHPLT ) 

DO 800 J-1,NPLT 
INET-0 

GO TO (300,301 *302, 303# 304) J 

300 IF ( IFLAGX.LT. 1) GO TO 299 
INET-1 

DO 38 I*l,NPTl 

IF(X( I ).GT.XHIN.AND. Yl< I ) .GT .YPIIN ) CALL TPLOT(X( I),V1(I),0,0) 
IF(X(I).GT.XPIIN.AND«Y1(I).GT. YPIIN) L-I 
IF(X(I ).GT.XWIN.AND.Y1(I ).GT .YNIN) GO TO 299 

38 CONTINUE 

299 IF(INET.LT.l) CALL TPLOT(X( 1 ), Y1 ( 1 ),0,0 ) 

DO 39 I*L,NPT1 

IF(X( I ) .GT.XflAX.OR.Yl(I).GT.YPlAX) GO TO 39 
IF(X( I J.LT.XNIN.OR. Yl( I ).LT. YPIIN )GO TO 39 
CALL TPLOT(X(I),Y1(I),I,0) 

IFd.EQ.NPTl )CALL TPLOT(X(NPTl ), Y1 (NPT1 ),0,0 ) 

39 CONTINUE 
CALL VECTOR 
CALL POINT 

CALL TPLOT(0. 0,0. 0,0,0) 

GO TO 800 

301 IF ( IFLAGX.LT. 1) GO TO 400 
INET*1 

DO 402 I*1,NPT2 

IF(X(I).GT . XWIN. AND. Y2( I ) .GT. YMIN ) CALL TPLOT(X( I),Y2(I),0,0) 
IF(X(I ) .GT . XFIIN .AND. Y2( I ) .GT . YNIN) L*I 
IF(X( I ) .GT.XNIN.AND. Y2( I ) «GT. YPIIN ) GO TO 400 
402 CONTINUE 

400 IFCINET.LT. 1) CALL TPLOT(X( I ),Y2( 1 ),0,0) 

DO 40 I«L,NPT2 
CALL VECTOR 
CALL CH0UTC2?) 

CALL CH0UTC97 ) 

IF(X(I).GT. XNAX . OR . Y2 ( I ) . GT . YN AX ) GO TO 40 
IFCXCI ).LT.XHIN.0R.Y2(I ).LT.YNIN) GO TO 40 
CALL TPLOTCXCI ),Y2(X ),I,0) 


IFCI.EQ.NPT2) CALL TPLOTCXCNPT2),Y2CNPT2>,0*0) 

40 CONTINUE 
CALL POINT 

CALL TPLOTC0.0 # 0.0,0 # 0) 

GO TO 800 

302 IF( IFLAGX . LT « 1 ) GO TO 404 
INET-1 

DO 403 I-i,NPT3 

IF(XCI ).GT.XHIN. AND.Y3(I KGT.YHIN) CALL TPLOTCXCI ) # Y3C I ) # 0,0) 
IFCXCI ).GT.XHIN. AND. Y3C I >.GT. YHIN) L-I 
IFCXC I )«GT .XHIN. AND.Y3C I ).GT. YHIN) GO TO 404 

403 CONTINUE 

404 IF ( INET * LT « 1 ) CALL TPLOT(X< 1 ),Y3C1 >,0,0) 

DO 41 I-L,NPT3 

CALL VECTOR 
CALL CH0UT(27) 
rail rumiTf qqi 

IFCXCI ).GT .XHAX.OR. Y3( I ) .GT . YNAX) GO TO 41 
IFCXCI). LT. XHIN. OR. Y3CI).LT. YHIN) GO TO 41 
CALL TPLOT CXCI)»Y3CI) # I,0) 

IFCI.EQ.NPT3) CALL TPLOTCXCNPT3),V3CNPT3),0.O) 

41 CONTINUE 
CALL POINT 

CALL TPLOTC6.0 # 0.6 # 0 # 6 ) 

GO TO 800 

303 IF ( IFLAGX.LT. 1) GO TO 406 
INET-1 

DO 405 I-1.NPT4 

IFCXCI ). GT.XHIN. AND. Y4CI ).GT.YHIN) CALL TPLOTCXCI ),Y4CI ) # 0,0> 
IFCXC D.GT.XHIN. AND. Y4CIJ.GT. YNIN) L-I 
IFCXCI). GT.XHIN. AND.Y4CI ).GT. YNIN) GO TO 406 

405 CONTINUE 

406 IFCINET.LT.l ) CALL TPLOTCXCI ).Y4C1 ).0,0> 

DO 42 1-L.NPT4 

CALL VICTOR 
CALL vH0UTC2?) 

CALL CHOUTC100) 

IFCXCI ).GT. XHAX.OR. Y4C I ).GT. YNAX) GO TO 42 


IF(XtI).LT.XI!IN.0R.Y4(I>.LT,YHIN) GO TO 42 
CALL TPLOTCXCI ) # Y4CI ) # I # 0) 



72 


Cttt DATA PLOTTING IS COMPLETE, LABLE SCALES ON AXES 
CALL ALPHA 

call scALE(i.e # i.e # e.e # e.e) 

CALL CHSIZEO) 

XO-XOS 


100 


66 

961 


► OKU AT C 1 H* # tt 1 J . 6 1 

YL«YL+TICY*YFACT 

STY-STY+TICY 

IFCSTY.LE.TY) 60 TO 56 

XT-STX+RLX 

VL»YOS-20.0 

XL-XOS-60.0 

XLZ-XL-120. 

CALL TPLOT (XL, YL, 0, 0 ) 
IF((XL-XLZ).GT.11S.0) URITEt IS, 901 ) STX 
F0RHAT(1H*,G13.2) 

XL-XL+TICXHXFACT 

STX*STX+TICX 

IF(STX.LE.XT) GO TO 60 

CALL UECTOR 

RETURN 

END 


PIP> 


SUBROUTINE UATNXT(IR) 

CALL ERASE 

xs-se. 

YS-500. 

CALL TPLOT(XS,YS,O,0 ) 

WRITE! 1,20 ) 

20 FORHAT(lH+,'DO YOU UISH TO*') 

YS-YS-25 • 

CALL TPLOT(XS,YS,0,0 ) 

WRITE(1,23) 

23 FORMAT (1H*,'1 HAKE OTHER PLOTS FROH SAHE 

* FILE WITH SAHE TITLES' ) 

YS-YS-25. 

CALL TPLOT(XS,YS,0,0) 

WRITE (1,24) 

24 F0RHAT(1H*,'2 CALL SCALING OPTIONS') 
YS-YS-25. 

CALL TPLOT(XS,YS,0,0) 

WRITE (1,25) 

25 F0RHAT(1H+, '3 RETRIEVE ANOTHER DATA FILE') 

YS-YS-25. 

CALL TPLOT(XS, YS,0,0 ) 

WRITE (1,30) 

30 F0RHATdH*,'4 HAKE OTHER PLOTS FROH SAHE 
t FILE WITH NEW TITLES') 

YS-YS-25. 

CALL TPLOT(XS,YS,0,0) 

WRITE (1,35) 

35 FORHATdH*, '5 STOP') 

YS-YS-25. 

CALL TPLOT(XS.VS,0,0> 

READ( 1,22 ) IR 
22 F0RHATCZ2) 

RETURN 

END 


PIP> 






SUBROUTINE XLABEL (XS,VS, ARRAY, NPA ) 

Cttt THIS ROUTINE OUTPUTS TEXT IN THE *X DIRECTION ON THE 4014 SCOPE. 
CXtS XS&VS ARE ABSOLUTE SCREEN COORDINATE LOCATIONS IN TEKPOINTS, 

C<» ARRAY IS THE FIRST LOCATION IN WHICH THE TEXT ARRAY IS STORED 

C*« (4 CHARACTERS/LOCATION), NPA IS THE NUHBER OF CHARACTERS TO BE 

Cttt OUTPUT TO SCOPE LUN-IS. 

COHHON/LUN/IS, IP, IC, ISD 
DIMENSION ARRAY(l) 

CALL SCALEtl. 0,1. 0,0. 0,0.0) 

CALL ALPHA 

CALL TPLOTCXS VS 0 0) 

cm NOTE: PDP ISSUES* LINE FEED COMMAND BEFORE WRITING, USER MUST 

cm ADJUST YS ACCORDINGLY, THUS FORMAT CONTAINS 1H* AND VS IS START 

NW*NPA/4 

IF((NPA-NW*4).NE.0) NU-NW+1 
WRITE! IS, 100) ( ARRAY (I),I a l,NW) 

100 FORMAT C 1H+, 30 A4) 

RETURN 

END 


e» 


PIP> 


SUBROUTINE VLABELCXS, VS, ARRAY, NPA ) 

Cttt THIS ROUTINE UILL UERTICALLV PRINT A CHARACTER STRING ON THE 
Cttt T4014 STARTING AT ABSOLUTE LOCATION (XS,YS) IN TEICPOINTS. DATA 

Cttt TO BE PRINTED ARE STORED IN ARRAY UITH 4 CHARACTERS TO THE WORD. 

CXSt NPA IS THE NUilBER OF CHARACTERS TO BE PRINTED. 

CONHON/LUN/IS, IP, IC, ISD 
DIHENSION ARRAY ( 1 ),X( 1 ) 

LOGICALX1 IX(4 ), IBLK 
EQUIVALENCE (X(1),IX(1)) 

DATA IBLIC/IH / 

NU-NPA/4 

IF ( ( NUT4-NPA ) . EQ . 0 )GO TO SO 
X(1 >• ARRAY (NU+l) 

DO 10 1-1,4 

IF( CNUX4+I ) .LE.NPA) GO TO 10 
IX(I)-IBLIC 
10 CONTINUE 

ARRAY(NU+1)-X(1) 

NU-NU+1 

SO CALL SCALEU. 0,1. 0,0. 0,0.0) 

CALL ALPHA 

CALL TPLOTtXS, YS, 0, 0 ) 

IC-1 

30 XCl)-ARRAY(IC) 

DO 40 1-1,4 
CALL CH0UT(IX(D) 

CALL CHOUT(IO) 

CALL CH0UTC8) 

40 CONTINUE 

IC-ICM 

IF( IC.LE.NU )G0 TO 30 

RETURN 

END 


PIP> 


CHPLOT VARIABLES 


ANPLT- 

CRLX- 

CRLY- 

CRLYR- 

IC- 

IFILE- 

IFLAGX- 

IFLAGY- 

IG- 

IGX- 

INET- 

INRLX- 

INRLY- 

INRLYR- 

IR- 

IS- 

ISTOPF- 

ISTOPT- 

IVECT- 

K- 

L- 

NCH(1. . 
NPLT- 
NPT- 
NPT(1. . 


number of plots on the graph (a real number) 

Integer length of X axis 
Integer length of Y axis 

variable used in checking plotting efficiency 
loop counter for the number of words to be printed 
plot file name (29 characters) 

determines whether to reset limiting X values (O-yes , 1-no) 
determines whether to reset limiting Y values (O-yes ( 1-no) 
grid line variable (-1 for grid lines, 1 for no grid lines) 
indicator for grid annotation ( .LT. 0-grid lines, .GT. O-tic marks) 
determines if limiting values have been reset (1-yes, 0-no) 
integer number of divisions on the X axis 
integer number of divisions on the Y axis 
variable used in checking plotting efficiency 
represents the next plotting option 

used as the unit number that data is read from or written to 

new file decision (0-yes, 1-no) 

new title decision (0-yes, 1-no) 

determines whether to call VECTOR (0-no, .GT. O-yes) 

counter for which channel is being plotted 

counter for point being plotted 

- channel numbers to go on the plot 
number of plots on the graph (max-5) 
number of points on the plot 

- number of points on the plot for five channels 
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NRDAT- associate variable for opening file 

NSCH- number of channels at the site 

NT- number of channels at the site + 1 for time 

NW- the number of words to be printed (4 characters to a word) 

RLX- length of the X axis 

RLY- length of the Y axis 

RLYR- variable used in checking plotting efficiency 
RS- variable for screen location 

STX- starting X value on axis 

STY- starting Y value on axis 

STYR- variable used in checking plotting efficiency 
T1CX- length between tic marks on the X axis 

TICY- length between tic marks on the Y axis 

TICYR- variable used in checking plotting efficiency 
TPX- number of tekpoints available in the X direction 

TPY- number of tekpoints available in the Y direction 

TY- final value on the Y axis 

XFACT- X scaling factor (tekpolnts/scaled unit) 

XL- variable used for plotting (used in calling TPLOT) 

XLZ- used in labeling axes after data plotting is complete 

XMAX- maximum abscissa value 

XM1N- minimum abscissa value 

X0- origin abscissa location for the grid in absolute screen coordinates 

referenced from screen (0,0) 

X0S- used in checking plotting efficiency (used in calling SCALE) 

XS- T4014 abscissa location 

XT- final X value on axis 


78 


YFACT- Y scaling factor {tekpoints /scaled units) 

YFACTR- used In checking plotting efficiency 
YL- used for plotting (used in calling TPLOT) 

YMAX- maximum ordinate value 

YMIN- minimum ordinate value 

Y0- origin ordinate location for the grid In absolute screen 

coordinates from screen (0,0) 

YOM- used in checking plotting efficiency 

YOR- used in checking plotting efficiency 

YOS- used to check plotting efficiency (used in calling SCALE) 
YS- T4014 ordinate location 


CHPLOT ARRAYS 


DATA(IOO)- contains data for one record 

IFILE(30)- plot file name 

IFILET(50)- plot caption 

PID(4)- plot caption for left side 

SAV(1. .5) (300)- contains channel data for one day 

SX(1. .5) (300)- units for the channels on the plot 

TIME(300)- stores time data for one day 

X(3Q0)- contains X values of points for one day 

XTITL(4)- contains the units for the independent parameter 

Y(l. .5) (300)- contains Y values of the points for 5 channels for one day 

UNITS(4)- contains the units for the dependent parameter 
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APPENDIX E 


WRTCNL SOURCE LISTING 


blank 




t!()T FKltf® 


SUBROUTINE INPUTC ICHNLS,NCH,NOUT,XNAMES,XNUN ) 

DIMENSION ICHNLS(1 ),XNAMES(1 ),XNUM(1 ) 

LOG I CALX 1 IFILEC30) 

CONHON/FILE/IFILE 
WRITE (1, 10) 

10 FORMAT (T2, 'ENTER NO. OF CHANNELS AT SITE ' > 

READ(1,15) NCH 
IS F0RNATCI2) 

WRITE (1,20) 

20 FORMAT (T2, 'ENTER THE NO. OF CHANNELS TO BE OUTPUT (10 MAX)') 

READ(1,15) NOUT 
WRITE (1,25) 

25 F0RMAT(T2, 'ENTER THOSE CHANNELS NOS. AS: 01,02,03, ') 

READ(1,30) (ICHNLS(I),I-1,N0UT) 

30 FORMAT (10(12, IX)) 

WRITE (1,35) 

35 F0RHATCT2, 'ENTER MEASUREMENT TYPE , AND NUMBER, AS: 

rmoe,EPiei,uuioo, '/' note: program expects 2 letter 

„ * IDENTIFICATION OF MEASUREMENT TYPE-'/' SINGLE LETTER MEASUREMENT 

® t TYPES SUCH AS TEMPERATURE'/' SHOULD BE INPUT AS TT,UW,ETC.') 

READ (1,40) ( XNUM (I ), XNAHES ( I ) , I - 1 , NOUT ) 

40 F0RHAT(10(A2,A3,1X) ) 

URITE(1,4S) 

45 FORMAT (T2, 'ENTER THE DATA FILE NAME OF THE STORED DATA') 

READ(1,50) (IFILE(I), 1-1,29) 

50 F0RMATC29A1) 

RETURN 

END 


PIP> 


SUBROUTINE OUTPUT ( ICHNIS,DATA, XNAHES* XNUR, NOUT, NCH ) 
OINENSION ICHNLSd ),XNANESd ),XNUHd ) # DATA(1 )«DAT(96) 
LOGICALtl IFILEOO) 

COnnON/FILE/IFILE 

J-l 

NRDAT-1 

ICNT-1 

URITE(3,3) IFILE 
3 F0RHAT(5X,29A1 ) 

URITEO, 1 ) (XNUHd ) # XNANESd ), 1*1, NOUT) 

1 F0RHAT(5X, 'TIHE'«4X, 10(3X,A2,A3,3X)//) 

2 READ ( 2 ' NRDAT , END *30 , ERR * 30 ) (DAT(I ),I-I,NCH*1 ) 

DO 10 1*1, NOUT 

Ji*ICHNtSd)M 
10 DATA(I)-DAT(JJ) 

IFCICNT .LE. 50) GO TO 20 
URITE<3,4> 

4 F0RNAT(1H1 ) 

URITE(3,3) IFILE 

URITE(3,1 ) CXNUHC I ),XNAHES( I ), 1*1, NOUT ) 

ICNT-1 

20 NRDAT-NRDAT+1 

WRITE* 3, 15) DATdMDATAd). 1*1, NOUT) 

15 FORHAT(F10.4,10F11 .2) 

ICNT-ICNTM 

J*J*1 

IF( J.GT.270) GO TO 30 
GO TO 2 

30 CL0SE<UNIT-3,DISP0SE- 'PRINT ',ERR-900) 

900 RETURN 
END 


PIP> 


VRTCNL: 


909 

PIP> 


DIMENSION DATA ( 10 ) # ICHNLS ( 10 ) , XNANES ( 10 ),XNUN( 10 ) 

LOGICAL! 1 IFILEC30) 

CONHON/FILE/IFILE 

OPEN C UNIT- 3. NAHE - ' PRT . LST ' , TYPE ■ ' NEW ' , ERR-900 ) 

CftLi ERASE 

CALL INPUT ( ICHNLS ( 1 ) , NCH, NOUT, XNANES Cl), XNUN Cl)) 

NT-NCH+1 

0PEN(UNIT-2,NANE- IFILE, TYPE- 'OLD' , ACCESS- 'DIRECT' , 
*RECORDSIZE-NT,ASSOCIATEVARIABLE-NRDAT, ERR-900) 

CALL OUTPUT ( ICHNLS Cl), DATA Cl) # XNANES ( 1 ) , XNUN ( 1 ) , NOUT # NCH ) 
CLOSECUNIT-2, ERR-900 ) 

STOP 

END 


WRTCNL VARIABLES 


ICNT- 

NCH- 

NOUT- 

NRDAT- 

NT- 


line counter (50 lines/page) 
number of channels at the site 
number of channels to be output 
associate variable for opening IFILE 
number of channels at the site + 1 for time 


WRTCNL ARRAYS 


DAT(96)- 

DATA(l)- 

ICHNLS(IO)- 

IFILE(30)- 

XNAMES ( 10 ) - 

XNUM(10)- 


array to store data for time and channels 

stores data tor each channel to be output 

array of channel numbers to be output 

data file name of stored data 

array of measurement number 

array of measurement type (TT,WW,etc.) 
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